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NOMENCLATURE 

f(t)  any  function  of  time 

Zf(t)  a  train  of  impulses  (equally  spaced,  T  seconds 

apart)  whose  areas  are  modulated  by  f  (t)  at  the 
sampling  instants 

f (s)  the  Laplace  transform  of  f (t) 

Z?(s)  the  Laplace  transform  of  2f  (t) 

,.  the  largest  absolute  value  of  error  in  the 

fixed  period  of  computation  (0,8  seconds) 

T  sampling  interval  size 

Error  (exact  solution  -  approximate  solution) 

Relative  Error  (  Error  /Exact  solution)  multiplied  by  100.0 

I  the  n*"  terra  of  an  ordered  sequence 

u  "approximate  equality* 


nrraoaicTioH 

The  purpose  of  thin  report  is  to  discuss  Mat  aspects  of  the 
approximate  solution  of  differential  equations  using  Newton-Cotes 
quadratures,     The  goal  is  to  illustrate  a  relatively  straightforward 
method  for  approxir.^.tiiic  U.e  solution  of  both  linear  and  non-linear 
differential  equation',,.      .11 3  aethod  is  suitable  for  direct  application 
on  a  digital  computer, 

2^transfemattto  la  Introduced  in  the  solution  of  differential 
equations  as  a  "bookkeeping"  procedure  for  keeping  track  of  desired 
coefficients*    Kalijak'a  "Irapezoidal-Gonvouition"  (1)  and  Simpson's 
Rule  are  the  — nhawlani  for  obtaining  the  approximate  analysis  of 
differential  equations. 

THE  ROLE  OF  2UTR&NSFORM3 

If  ideal  samplers  are  admitted,  then  ^transformation  pemits  an 
organised  approximate  solution  of  a  linear  system  at  the  sampling 
instants.    In  the  case  of  a  non-linear  system  the  solution  cannot  be 
effected  solely  by  the  algebraic  manipulation  of  ^transforms,  but  the 
solution  can  be  implemented  by  non-linear  recurrence  relations.     The 
solutions  thus  obtained  are  always  approxiaations  of  the  true  value  at 
the  sampling  Instants. 

If  an  arbitrary  function  of  time  f  (t)  is  always  multiplied  by 
g(t)  defined  by 

g(t)  =  J  oT(t-nT)  (1) 


r*»0 


where  T  is  the  time  duration  between  two  successive  pulses,  then  tha 
resulting  product  can  be  expressed  in  Laplace  transforms  as 


.C<t)  f(t)l  •*!*><$  fU)i 

L      J   2TTJ  J9 


■ahoy© 


&•>  -  L 


£>£ 


oT(t-nT) 


f(A)g(s-A)dA 


tW^**«»  ...  4e~»*T 


1 


(2) 


(3) 


|e*>* 

and  R  is  a  region  in  the  complex  A  -plane  containing  the  poles  of  f (A). 

SMMWMPt 


t<t)  f(t)l  •  JL'A  f(x) — 4 

L  J        2TrjX  1-e-T( 


3- A) 


dA 


w 


This  particular  Laplace  transform  is  the  transform  of  the  sampled  function 
f(t). 

Xqtttlon  (4)  can  then  be  written 


df 


:ttj 


l/(1-e*r**"  ')  can  be  developed  as  a  geometric  progression  with  ratio 
0~T(s-X)  ^  ^  ^tiai  term  of  one. 
Equation  (5)  is  then  written 


[s»A) 


ft2TTj      JT 

SS27TJ  / 


(7) 


The  integral  represents  the  inverse  Laplace  transform  and  we  obtain 


Zf  -   >  e"nTs  f  (nT) 


(8) 


n=0 


,-Ts 


We  then  define  e~iS  =  s,  f (nT)  s  f  and  equation  (8)  can  be  shortened  to 


*  ■  2  n; 

n-0 


(9) 


NUMERICAL  APPROXIMATIONS  OF  THE  CONVOLUTION  INTEGRAL 

The  solution  of  a  linear  differential  equation  with  constant 
coefficients  can  be  considered  as  a  convolution  problem. 


Left  Riemann  Sum  Approximation 


f(r)g(t~?)dT=J 


f(r)g(t-7)dT=l      g(T)f(t-T)dT 


is  equal  to 

H)T  n-1      r(K+1)T 

f(T)g(nT-T)dT=y     / 


i-1      f( 

2/ 


k=0     KT 

and  is  approximated  by 


2/ 

k=0     KT 


(10) 


:-i 


f k6n-k  "  T  L  gkfn-k 

k=u 


The  above  becomes 

n-1  ■ 


T  2  '#**  •  T  2  f  kgn-k  "  ^0 
k=0 


(11) 


(12) 


Z- transformation  yields 
n 


2  f ^n-k  "  f  ne 


=  T(Zf)(Zg)  -  TZfg0  =  I 


f| 


(13) 


If  we  use  the  symmetric  property  of  the  convolution  integral  and  take  the 
average  of  the  two  results  we  find  that 


f* 


T(Zf)(Zg)  -|(f0Zg  +gQZf) 


(1*) 


Right  Riemann  Sum  Approximation 


Ihe  right  Riemann  sum  employs  the  right  ordinate  and  yields 

,  n-1 

,  fk+1   ^n-k-1  =  T  2  gk+1  fn-k-1  (15) 


fk+1  KM  "  T  2  fk  ^k  '  T  2fk  ^k  -  Vn        (16) 


k=1 


k=0 


Z- transforming  (16)  we  have 


Z  T 


k=1 


fk  ^n-k  "  f0gn 


=  T(Zf)(Zg)  -  TfQzI 


07) 


If  we  use  the  symmetric  property  of  the  convolution  integral  and  take  the 
average  of  the  too  results  we  find  that 
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■  T(Zf)(Zg)  -1  (fQZg  +gQZf) 


(18) 


Trapezoidal  Convolution 


By  definition,  the  area  under  the  curve  using  trapezoids  is 


Jn  =  2  Z  f k  Sn-k  +  f 


n-k   *k+1  Bn-k-1 


(19) 


■ 


~1 


(20) 


k»0 
^transformation  yields 


2  2  *  ^  "  f»«0  +  2  2  fk  Sn-k  *  V» 


(21) 


:=0 


zjffej  - 1  (zf)(z£)-fz?  +  f<zf)(z£)  -  |r0zj 


nd 


Z  [fj]  »  T(£)(z£)  -  1(2^  *  Zjf0) 


(22) 


It  Is  observed  that  equations  (W),  (13),  and  (22)  are  all  the  same, 
by  using  the  qyaastrlo  property  of  the  convolution  integral  the 
equation  of  trapezoidal  convolution,  (22),  may  be  arrived  at  by  using 
loft  or  right  hand  Rieaann  sums. 


-  odifiod 


*s  Convolution 


^> 


y0 


r\ 


*.   k2 


-TOT  x 

Fig,  1.     Curve  used  for  Simpson's  rule. 

In  Fig,  1  we  know  the  ares,  between  the  ordinates  yQ  and  y2  is, 

T2  -  ft*  ♦  %  *  y2)  (23) 


ol  is  the  area  between  the  1th  ordinate  and  the  yQ  ordinate, 
T1<5  is  the  area  between  the  1th  and  3th  ordinates# 


rule  depends  on  an  oven  number  of  intervals;  this  is  one 
aspect  of  Simpson's  rule  whiah  oust  be  modified*  It  will  be  required 
that  only  one  interval  be  usod  in  an  approximation.  Therefore  wo  prooeed 
as  follows*  Simpson's  rule  requires  three  points  to  fit  the  parabola 
aX2  ♦  bX  ♦  o.  Jron  Pig.  1 


/*0  ° 

X  -  /     ttifantitt  »  *&  ♦  ^£~  ♦  «  &L  „  £2~  +  0i    (2*0 

1    •J-T  3  2  -T        3  2 


•tX»-T  y0  »  aT2  -  bT  +  o 

at  X  «  0  y,  *  c 

lkZ»T  y2  •  aT2  ♦  bT  +  o 

Dm  next  step  is  to  find  sons  combination  of  y„,  y^  and  yg  to  yield    0^. 

The  first  observation  is  to  multiply  this  partioultr  aombiaation  by  T  to 

obtain  the  desired  aT3,  bT2,  andoT.     Therefor,  equating    OJ  and  the  sua 

of  y0.  y1 .  and  y2  we  have 


T(By0  •  Gg1  +  Dgg)  ■  *t-  -  ^  ♦  OT 


(25) 


B,  C,  and  D  represent  the  desired  co«ffioients. 
Therefore,  we  have 

TCBaT2  -  BbT  ♦  Do  ♦  Co  ♦  DaT2  +  DOT  +  Do)  ■  &-  -  ^  +  eT     (og) 

3    2 

How  equating  ooofficionts  of  aT3,  bT2,  eT; 

(B+D)  -  1/3,  (-3+D)  ■  -1/2,  B+C+D  ■  1 

Solving  these  equations  we  find  thai;  I  *  5/12,  D  »  -1/12  and  C  ■  2/3. 

From  this 

X   S  <lV0  +  fc  "  ^2>T  (27) 


and 


we  find  04o:  We  know  that 

1 

Y12  -V2  .^  -  ^y0  ♦  *»,  ♦  J»,  .  ^0  -  £*,  .  -Uy2 

^2  "  <-  fo  +  ^1  +  ^2>T  <28) 

If  we  calculate  02o  tne  same  way  as  0. ,  we  get 

X        •  (.  ±y   +  %   +  JLjf  (29) 

•~->     12  ■   3  ■   II  * 

Also, 

X^  m  fr  J-u  ♦  %  ♦  -ly  )T  (30) 

**      12     3  J   12  * 

From  (28),  (29),  and  (30)  it  is  observed  that  a  pattern  exists  resulting 
in  the  following  recurrence  relation 

^(n-1 ) ,n  =  ("  ^n-2  +  Fn-1  +  ^»>T  ^  } 

n  =  2,  3.  ^.  5.  6,  ... 
From  (31)  we  can  establish  the  relation 

n  =  2,  3,  ^  5.  6,  ... 

The  preceding  development  now  enables  us  to  make  an  approximation  to 

the  convolution  integral, 

ZnT 

I   f(nT-T)  g(T)dT  (33) 

Jo 


(33) »  f(nT-T)  g(r),  represents  the  ordinate  values  of  yQ,  y^»  y2, 
7y   ...  eto.  for  various  values  of  n.  If  n=1  we  have  from  equation  (33) 

t 


f(X-T)  g(r)dT 


T  s  0{         *  •  fl8o 

f  »  2;  j2  «  f^ 

IJornaHy  f_-j  would  bo  taken  to  bo  ~ero,  but  In  this  ease  the  extra  point 
is  needed  tc  doteraine  the  direction  of  tha  curve,  since  S2jtp9on*s  rule  is 
based  on  a  ourve  passing  through  3  points. 
Next  it  is  known  that 


When  n=2  we  have  frora  equation  (33) 

t 


f(2T-T)  g(T)dT 


?*  1*  yt  »  f1g1 

f  •  H  y2  -  f0gg 

and  we  imov 

\  B  ^<y0  ♦  ^|  *  y2) 

therefore 

T2  •  5<*2*>  +  *1*1   +  W  <35) 


When  n=3  we  have  from  equation  (33) 
/*3T 


f(3T-T)  g(T)dT 
'0 


I  =  0,  vQ  =  f3g0 

t  =  ii         r,  ■  v, 

f  =  2,  y2  =  flg2 

From  equation  (32), 

i        ^  12  1       3  2      12  J 

_  1/ 


3(y0  +  ^1  +  y2}  +  (-  171  +  72  +  ^3)T 

'  f(f3S0  +  ^2*1  +  f1*2>  +  (-  Jfzh  +  fflS2  +  ^3)T     (36) 


When  n=4  we  have  from  equation  (33) 

At 

Jo 


f(^T-r)  g(T)dT 
'o 


T  = 0{         y0  =  Vo 

f»*l  y1  =  f3g1 

T  =  2j  v    =  f„g„ 

'2         2&2 

I  ■  3J  73  =  Iff, 

T  =  *  hm  foh 

From  equation  (32) 
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\  "  T^S0  +  fr*   +  7^2*0  +  ^1^3  +  l^W  (37) 

In  a  like  manner, 

Y5  "  T(^5S0  +  K%  +  T?3g2  +  f2*3  +  T?1g1  +  ll^  (38> 

and 

X6  =  T^6S0  +  ^5g1  +  W*M  +  f3g3  +  f2%  +  1^1*5  +  llfog6) 

(39) 
At  this  point  it  is  observed  that  a  definite  trend  is  being  followed 
in  that  the  first  three  co off icients  and  the  last  two  coefficients  of 
''if.*  Oct  0g»   •••  0„.  are  remaining  constant.     This  is  the  fact  that  allows 
a  closed  form  to  the  approximation  of  the  convolution  integral. 
Nejct,  write  each  area  in  the  form 

r0  =  o 

Y3 =  Wo +  &* +  "1% +  ^0^3  m 

\  ■  %*l*0  +  ^1  +  7^2%  +  ^163  +  T^OS*  (*3) 

^5  =  ^Vo  +  ^1  +  ^f3g2  +  Tf2*3  +  ^f1%  +  ^0S5  m 

\  *  ^0  +  ft 5*1  +  jf^2  ♦  Tf3g3  ♦  Tf^  ♦  l|Tflg5  ♦  ±Tf0g6 

m 


11 


It  should  bo  noted  at  this  time  that  tho  subscripts  of  each  ordinate 
value  add  arithmetically  to  n,  i.e.,  from  (41),  f^,  fjgj,  f^.  we  have 
each  ordinate  subscript  adding  to  2  which  Is  the  subscript  of  the 
corresponding X  Zt  is  also  observed  that  each  0^  has  n+1  terns  and  that 
all  terms  added  for  each  new^,  (for\,  0^,  0^,  »..)•  besides  the  first 
three  and  last  two  teats,  have  a  oonstant  ooof ficiont  of  1 •  One  more 
item  that  should  be  taken  into  account  is  the  fact  that  in  each  0^, 
(n  »  2,  3,  4,  5,  6,  .,,),  the  following  eequenoe  of  ordinates  occur i  for 
\,  we  have,  f^gn,  f^,  f.^,  f^,  f^;  Mai  the  subscript  on  f 
starts  at  n  and  decreases  to  f  0  in  n+1  terms  and  the  subscripts  on  g 
do  just  the  opposite, 

(ZfXzJ)  -  J  fk  ^ 

"  Voa°  +  (f1S0  *  Vl)a  +  (f2S0  *  W  W*  +  - 

(46) 

flBtmlim  (45)  with  (C  a  0  and  equations  (40)  through  (45)  we  see  that 
0q  corresponds  to  the  s  term  of  (Zf)(2g),  ^  corresponds  to  the  z  term, 
Ol  to  the  er  term,  etc. 

In  this  next,  step  ?  *Jill  be  suppressed  from  equations  (40)  through 

(45)  and  the  necessary  corrections  to  (Zf)(2g)  will  be  made  to  account  for 

coefficients  other  than  1 ,  or  any  extra  or  deleted  terms  that  need  taken 

caro  of  in  equations  (40)  through  (45), 

Comparing  (46)  with  equations  (40)  through  (45)  the  following 

corrections  are  applied  to  (2?)(^)| 

(2f)(2g")  -  ^,2?  ♦  Jg^Z?  -  ^z2^  ♦  ^aZg  -  -Z^Zg 


Collecting  toms  we  havot 

<^)(2S)  -  (§80  -  |fta  *  ^s2a2)zi  -  (-  Aqa  ♦  ir0)5  (<*?> 

Equation  (*f?)  is  a  corrected  equation  for  the  first  throe  terns 
Ust  two  terns  mentioned  on  page  11. 

Multiplying  the  last  two  terns  of  (ty?)  out 

f0   f1   f2   f3   *lt       f5 
§80   -fa    +  fe> 

§*o«o      fc%       ^2S0       fcfl,       ^gQ       §*5«o 

*Wft   -£figi   "Kgi   *^3S1   -K«i   ~i?fi 

l?0h      j?&      feg2      1^3S2      fo%      ^5g2 

and 

g0        «l        g2        63        g*        «5 

■ft       "12^1 

Ifai    #Vl         ll*0g2        ^0g3        #0%        1^0S5 

jfeS    "ffih    *jft%    «^1S3    -1^1%    -^1*5 


Adding  oooffioients  of  these  ttro  terns  and  inserting  *  minus  sign  wo  have 


13 


-  JiVo  -  Kei  -  3*2*0   -  ¥?0  ~  §Vo  -  1*5*0 

"  12*1*0      3*1*1      Kg1       5^1       JVl 
•  3*052  "  12*0*3"  12*2*2  *W  A 

"12*0%  "12*0*5 

JU  .     JU  . 

Wni  12^ 

Adding  these  coefficients  to  the  coefficients  of  (Zf)(zl)  we  get 

-  to*> +  (3*0*i +  to** +  <fc*>  *  Imi  +  Kg2)a2 

+  (3*3*0  +  &>*1   +  1^0g3)a3  +  (K*0  +  £f3S1+  Wjk  +  12*0% 
♦  fay  J*  ♦  C^A  •  j^   ♦  ^  ♦  ^  ♦  ||^  +  f2g3)a5 

<*a> 

Comparing  these  results  withX  •  0  and  equations  (40)  through  (4-5) »  the 
following  is  concluded! 

00  ■  0,  therefore  we  will  need  to  add  ^o^* 
For^,  we  need  to  add  -  -*-f  A»  a 

the  rest  of  the  0  •s  oorresporsd  exactly  and  no  further  correction  terms 
«jr«  required. 
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Therefore,  we  may  write  down  the  following  approximation  to  the 
convolution  integral. 

ftjS  *  T(zf)(zl)  -  KfgQ  -  JgiZ  +  ^2)zf  -  TC^q  -  ^flZ)zg 


&0h  "  12^-1  S2Z 


This  formula  will  be  called,   "The  Modified  Simpson's  Convolution  Formula". 


1 1  n 
1/s 


FOR  THE  INTEGRATOR  SUBSTITUTION  PROGRAM 


1/s11 


from  the 


Halijak  has  developed  an  approximate  Z- transform  of 
trapezoidal  convolution  formula  (1).     For  n  ■  1,  2,  3  this  approximation 
coincides  with  the  exact  solution. 

Approximation  of  Z|j/snJ  by  the  Modified  Simpson's 
Convolution  Formula 

In  a  similar  manner,  as  that  in  (1),  we  find  equation  (4-9)  yields 
zfi/sn+lj  mJ* 


Tz 


(1-z)2 


f*lp*ttlMY0* 


(50) 


n-r 


n  =  2,  3.  4.   ... 


IJhere 


rv  -    T    5-*6z-z2  ,  „        T    ,1     „    .    1,  „       1,  „2v 

°<  "  T2  ^TT"  •   »*  B  =  T2  (^1Z  +  T2g2Z  "  3V  ) 

Checking  for  accuracy,  equation  (50)  yields 


l/l 


„  J&m  m  .  tt.  +  i22  .  ^3  .  J^ 


2(1-z)J      36        18 


18 


* 


(5D 


as  compared  with  the  exact  solution 
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1/s3 


.   T2 


2(1-2p 


(52) 


This  rules  out  the  use  of  equation  (4-9)  as  a  method  of  approximating 


1/Sn 


,  for  it  should  at  least  be  exact  for  n=2  and  3  to  be  useful. 


The  Exact  Solution  of  Z 


l/sn 


Criswell  (2)  presented  an  algorithm  for  determining  the  Z-transforra 

of  tn/nl.  This  algorithm  generates  tables  of  Z-transforms  of  tn/nl. 

Table  1  of  Appendix  A  can  be  used  to  obtain  the  coefficients 

n-2 
A(n,p)  in  An(z)  =  ^  A(n,p)zp  for  the  exact  form, 


p=-1 


|  (1/sn)  =t^tn!^L,     n=1,  2,  3,  f,  ...  15 
(n-1)i  (1-z)n 


(53) 


Z-TRANSFORMS  USED  IN  THE  SOLUTION  OF 
DIFFERENTIAL  EQUATIONS 


Z-Transforms  Developed  with  Trapezoidal  Convolution 


The  following  equations  were  developed  by  Halijak  (1),  using 
equation  (22). 


i- 
a* 


2  1  1+2-2-       1  fl. 
~  2  1-z  8  "  2  1-z 


Z(l/sn)  = 


Tz 


T  1+z 


(Uzf[2  1-zJ 


n-2 


li 


1      T2; 


rr  1+z 


L   J    o-z)2b  1-z. 


n-2 


(Zg  -  jJBq) 


n  >  2 


n  >  2 


(5<0 
{55) 
(56) 
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Z-Transforms  Developed  from  the  Modified 
Simpson's  Convolution 


Using  Criswell's  exact  Z 


1/s" 


equation  (^9)  yields 


s  B 


1     . 

s7  S 


-  JL 

1-z 


ill + k  -  i^2)zi  -  iK  -  (ko  -  J* +  M  (57) 


(i-z)' 


(lis  .  V  +  J^3)Z|  . 
L12         6  12 


(3*o  ~  Tl"^)2  +  (5^1  "  P-z**2 


t* 


(53) 


Z  4r  I    *  -£-*  (U,  ♦  2,2  +  lz3       1    *)zg  .  (4,    +  -Lg  )z  - 

_s3  J     2(i-z)3l_i2      r      *       w  3       12 

(§g0  -  |ii  -  £g2)a2  +  (pi •  Wz3J 

In  general 

IL  |]  -  <*i7Z      I       +  1  _l      ,2<>  fjESt    An(2)  ) 

#    J  3*0       4«1Z       12  2         1(^)1      (1.z)n 


(59) 


(60) 


Z-Transforms  Developed  from  the  Averaged  Modified 
Simpson's  Convolution 


Using  the  symmetric  property  of  the  convolution  integral,  taking 
the  arithmetic  average  of  ths  results,  the  following  formulas  were 
developed: 


z 


17 
fi]  =  T(Z?)(z£)  -  f($r0  -  j^l  +  ^2z2)Zi  -  f^  -  3*1 8  #  lF2z2)^  + 

^oso-^g.i^+^i^)2  (61) 

1 1|  =  jl  [a  ♦  tt. .  -w  +  j*3)zg - a*  - 

Is    J       1-z    L8       24        24  24  8°° 

(J^o  -  Jet  +  ^s-1  +  fa^i  +  ^iz2j  (62) 

X  g]  .  -lL-   [A  .  ^2  +  1,3  .  1^  .  (^  +  ±g       .  J^)a  + 
js2    J       (1-z)2    I  6         12  3  12  3^       12  24  ■ 

4ft  +  Kl  -  fe^2  -  12s3]  (63) 

zri  5>-^ [A +  1*2*  «3.v  +  a,5)z5. 

Ls3    J      2(1-z)3  16        3  3         6 

<&o +  Ki +  2Wz  -  4%  -  £*i  -  ki  •  K)*2  + 

Since  there  were  no  original  equations  to  check  equations  (57)  through 
(64)  an  effort  was  made  to  find  some  simple  check.  It  was  observed,  after 
deriving  the  first  two  equations,  that  the  numerical  coefficients  of  Zg 
added  to  1  in  each  case  and  the  rest  of  the  numerical  coefficients  added 
to  -1/2.  Using  this  as  a  guide  it  was  found  that  this  was  also  true  for 
all  of  the  equations  developed  thereafter,  thus  enabling  a  simple  check 
for  accuracy. 

The  Averaged  Modified  Simpson's  Convolution  formulas  x^/ill  not  be  used 
in  the  numerical  analysis  since  they  contain  two  negative  subscripted 
variables,  of  which  only  one  can  be  found.  It  should  also  be  noted  that 
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formulas  eroato  recursive  relations  which  are  one  order  higher  than 
those  stained  by  the  Modified  Simpson's  Convolution  formulas*     This 
could  possibly  make  the  numerical  data  more  susceptible  to  round-off  error 
and  spurious  solutions* 

HUMERIGAL  EXAMPLES 

tte  following  differential  equations  are  solved  using  Trapezoidal 
Convolution  and  the  Modified  Simpson's  Convolution  formulas! 

j|  ♦  f  m  1  y(0)  «  0 

♦  ^t4y3l  y(0)-y(0)-o 

<&  +  2&L  +  &L  +  y  .  ,         ^(0)  »  y(0)  «  y(0)  «  o 
dp     dt2      dt 

dt     ^ 

Only  tho  equations,  the  type  of  solution  used,  the  reeurrenoe  relations, 
and  pertinent  data  are  shown.     The  step-by-step  solutions  and  computer 
progtMs  are  shown  in  the 


d 


Empirical  Error  Analysis 

A  comparison  of  the  computational  orrors  generated  by  various  •—^■g 
when  using  Trapezoidal  Convolution,  tho  ibdified  Simpson's  Convolution 
and  Simpson* s  1/3  rule  is  presented  ftp  the  first  order  differential  equation. 
Tho  data  for  Simpson's  1/3  rule  was  taken  from  Loi's  Master's  Report  (11). 
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va  Formulas 


a)     Solution  of  Sfc  +  y  »  1f  y(0)  ^  1     (Appendix  B) 
dt 

Sz&ot  solution, 

y(t)  *  1-e"*  (65) 

Trapeeoid&l  0>>ivo3.ution, 

yn  a . j-  -. 

Modified  Simpson's  Convolution,  -^ 

7n  a  — £ a — 

(1  +^-T) 


WT       {^}-<0#T<1   ♦&*  ^).T,T,  ...  I 


(67) 


T(1  +  J)  2T 

yi  a 3*  '      y2  "         ,1 

1-HP*  C  1"*T*  ^ 

3  3 

b)    Glutton  of  B  ♦  Ik  •  y  *  1 .  y(0)  a  y(0)  =  I    (Appendix  C) 


solution, 
y(t)  »  1-e~W"1  (68) 

Trapecoidal  Convolution, 

n      v        ^n-1       v       /Jn-k: 


y    _  _____ 

n  (TH) 


{%i}H°'  2'  1»   ••   1»   —f 


Modified  Siiapson's  Convolution, 

JL  -  (U£  -  2  ♦  |)y        -  (1   -  £  -  2t)v        -  T(X  +  iw 
^       ^  12                2yn-1       u        6       2i;yn-2       nl2      6;yn-3 
yn  . 

(1  +f)  (70) 

I  nj  2        12       26  4626 


where 


7}   • 

*<k- 

3       6; 

y    —                      - 

,+2I+f  r 

*                        h«2       ™A 

1   +  2T  +  —      ~ 
3         9 

c) 

Solution  of 

^  +  2&  +  2&  +  y  =  1ty(0)  =y(o)  =y(o)  = 

dt3      dt^        dt 

(Appendix  D) 

0 

Exact 

solution, 

y(t)  i 

=  l-e-t  -  te"1  - 

t2/2  e-t 

Trape 

soidal 

Convolution, 

Y-n 

n 

^3 

-  > +  f  V,  - 

o-f^-'f-1*. 

-3 

(71) 


(1  +  f )  r  (72) 

{\M0»  kl>  1'  1»  1»  ••• 

Modified  Simpson's  Convolution,  -> 

»  m  [x     -  (i22-  .1.3+  i22Lw       .  (3  .  3-r  .  152-  +  21-w 

r»      \y        V  4  2      '        24  ;yn-1       U      *        4  8  ;yn-2 

-  cf  - 1  +  ^  ♦  T>v3  +  4  +  f  +  fer**){  Ci  ♦  f )     (73) 
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I  ajl  4-         I         I         i  2^y,%  I  2  2 

-  2EL  .  &  «►  I&to  ^(T2  _  a  .  S3    &  x  y2  _  ^2  ^  *V2) 

4-3  8  4  <*  8  4         2  5 

T3,  T3,  T3.   ...I 
where  -^ 

d)  Solution  of  ^  +  y2  =  1 ,  y(0)  =0  (Appendix  E) 

Exact  solution, 

y(t)  »  1  -  -2—  (7^) 

1-K>2t 

Trapezoidal  Convolution, 


-1+^- 


y  = 3 n  ar  1     (j$) 

T 

Modified  Simpson's  Convolution, 

V         T* Tn-1       yn-1       12  V2         ; 


yn3 


where 

v  =  T ,  £  +  z£ .  5.67T7 

y1  3         15        105 

2T      8  3      64T£      226£ 
y2  31         15         105 
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Results 

The  error  measurement  for  the  first  order  equation  is  presented  in 
Table  2,  and  graphed  in  Figs.  2  through  6.  Decimal  place  accuracy  for 
various  sampling  intervals  for  the  second  order,  third  order,  and  non- 
linear differential  equations  is  presented  in  Tables  3  through  14. 


Tab; 

max     dt   » 

■  U   y(0)  =  0,  solution 

time  (t)  =  .8  sec. 

T 

Trap.  Conv. 

Mod.  Simp.  Conv. 

Simpson's  1/3  llxde 

.001 

.00000015 

.00002292 

— 

.002 

.00000025 

.00000496 

«• 

.003 

.00000110 

.00000092 

.00000038 

.004 

.00000055 

.00000567 

.00000039 

.005 

.00000070 

.00000202 

.00000014 

.006 

.00000087 

.00000050 

.00000027 

.007 

.00000097 

.00000325 

.00000031 

•  1w'v'^ 

.00000142 

.00000122 

.00000015 

•  009 

.00000235 

.00000031 

.00000014 

.010 

.00000253 

.00000224 

.00000012 

.012 

.00000389 

.00000029 

.00000014 

*m 

.00000582 

.00000064 

.00000008 

.016 

.00000738 

.00000129 

.00000018 

.013 

.00000960 

.00000031 

.00000015 

.020 

.00001191 

.00000036 

.00000015 

.022 

.00001442 

.00000084 

.00000015 

.024 

.00001708 

.00000034 

.00000013 

.026 

.00001996 

.00000013 

.00000014 

.030 

.00002332 

.00000046 

.00000015 

.00002654 

.00000051 

.00000014 

.032 

.00002269 

.00000040 

.00000011 

.034 

.00003443 

.00000023 

.00000011 

.036 

.00003863 

.00000074 

.00000011 

.038 

.00004320 

.00000053 

.00000011 

.040 

.00004779 

.00000049 

.00000010 

.042 

.00005278 

.00000110 

.00000013 

.044 

.00005778 

.00000094 

.00000009 

.046 

.00006300 

.00000092 

.00000013 

.046 

.00006341 

.00000154 

.0000001 6 

.050 

.00007475 

.00000154 

.00000006 

.052 

.00008059 

.00000150 

.00000012 

.054 

.00008624 

.00000213 

.00000006 

.056 

.00009346 

.00000216 

.00000008 

.058 

.00009941 

.00000220 

.00000008 

♦060 

.00010724 

.00000285 

.00000014 

.062 

.00011324 

.00000294 

,00000010 

.064 

.00012163 

.00000302 

.00000008 

.066 

.00013022 

.00000374 

.00000008 

.068 

.00013646 

.00000383 

.00000010 

.070 

.00014559 

.00000405 

.00000007 

.072 

.0001 5499 

.00000478 

.00000009 

.074 

.00016121 

.00000483 

.00000012 

.076 

.00017113 

.00000514 

.00000012 

.078 

.00018138 

.00000592 

.00000015 

Table  2.     (cont.) 


T 

Trap.   Conv. 

Mod.  Simp.  Conv. 

Simpson's  1/3  Rule 

.080 

.00019174 

.00000632 

.00000006 

.082 

.00019735 

.00000634 

.00000012 

.034 

.00020835 

.00000719 

.00000009 

.086 

.00022012 

.00000759 

.00000010 

.038 

.00023164 

.00000804 

.00000015 

.090 

.00023673 

.00000847 

.00000013 

.092 

.00024394 

.00000390 

.0000001 1 

.09^ 

.00026122 

.00000951 

.00000010 

.096 

.00027391 

.00001038 

.00000015 

.098 

.00023673 

.00001104 

.00000022 

0.100 

.00029981 

.00001169 

.00000016 
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10" 


10 


-7 


IV 


•Mod.  Simp.  Conv. 


I   \ 


Trap.  Convi 


.005 
T(sec.) 


.010 


Fig.  2.  Error  .curves  of  ™C  +  y  ■  1 ,  y(o)  =  0,  solution  time,  (t)  =  .8  sec. 
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/  > Trap.  Conv. 

10-5 

/ 

x~ 

f 

rt 

.s 

I 

« 

1 

1(T6 

/  i     n                  / 

/  •    /  «            i  \  ■  j 

1    /    !    /    »                i    \ 

1   /        /     i       r\      i     V 

■ 

1    /              i           l          '         *       |  \         Mod.   Simp.   Conv. 

•         u 

» 

t  1 
i  i 

u 

icr7 

-  >          i 

.01                .02                .03                 .04              .05 

T(sec) 

Pig.  3. 

Error  curves  of  §J  +  y  =  1,  y(0)  =  0,   solution  time,   (t) 
at 

=  ,8  sec. 
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10 


-6 


w 


10 


-? 


o-3— tf- 


.01 


Conv. 


Simpson's  1/3  Rule 


.02       .03 
T(sec.) 


.04- 


.05 


KLp-.  *i<.  Error  curves  of  -^+  y  =  1,  y(0)  =  0,  solution  time,  (t)  =  .8  sec. 

dt 


to"*- 


W 


10-5  _ 


10 


-6 
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Trap.  Conv. 


,'*S. 


Mod.  Simp.  Conv. 


^y 

^** 


/ 


* 


i    i 


0.05    o.o6 


0.07      0.08 


0.09 


0.10 


dy 
Fig.  5.  Error  curves  of  ,.  +  y  ■  1 ,  y(0)  =  0,  solution  time,  (t)  =  .8  sec. 
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Mod.  Simp.  Cor. v. 


■#—*■ 


0.05       0.06 


0.07 
.  T(sec.) 


0.08 


0.09       0.10 


Fig.  6.  Error  curves  of  ^  +  y  =  1 ,  y(0)  =  0,  solution  time,  (t)  =  .8  sec. 
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Table  3. 

'{ 2  )  Y/DT ( 2 )+2DY/DT  +  Y  =  l» MODIFIED  SIMPSONS  CONV. 


TIME 

EXACT 

APPROX. 

ERROR 

REL.  ERROR 

.003 

.00000439 

.00000449 

.00000010 

2.38955360 

.004 

.00000787 

.00000798 

.00000011 

1.41212190 

.005 

.00001244 

.00001246 

.00000002 

.16816080 

•  006 

.C0001789 

.00001793 

.00000004 

.20215180 

•  007 

.00002433 

.00002439 

.00000006 

.23626200 

.  ")08 

.00003174 

.00003183 

.00000009 

.27049206 

•  009 

.00004024 

.00004026 

.00000002 

.05330745 

•  oio 

.00004960 

.00004967 

.00000007 

.13499325 

.oil 

.00006004 

.00006006 

.00000002 

.03535487 

.012 

.00007134 

.00007143 

.00000009 

.12430682 

.113 

.00008371 

.00008377 

.00000007 

.07926767 

.014 

.00009703 

.00009709 

.00000006 

.05800132 

.015 

.00011132 

.00011138 

.00000006 

.05526410 

.016 

.00012656 

.00012664 

.00000008 

.06409507 

.017 

.00014286 

.00014287 

.00000002 

.01127701 

.018 

.00016000 

.00016007 

.00000007 

.04346^50 

.019 

.00017819 

.00017823 

.00000004 

.02215039 

.020 

.00019733 

.00019735 

.00000003 

.01416438 

.021 

.00021740 

.00021744 

.00000004 

.01703764 

•  022 

.00023841 

.00023848 

.00000007 

.02843805 

.023 

.00026046 

.00026048 

.00000002 

.00702598 

.024 

.00028344 

.00028343 

.00000000 

.00260724 

.025 

.00030735 

.00030734 

.00000001 

.00331867 

.026 

.00033219 

.00033220 

.00000001 

.00419644 

.027 

.00035805 

.00035801 

.00000004 

.01114379 

.028 

.00038472 

.00038476 

.00000004 

.00997341 

.029 

.00041242 

.00041246 

.00000004 

.00990010 

.030 

.00044113 

.00044110 

.00000003 

.00667146 

.031 

.00047065 

.00047069 

.00000003 

.00715600 

•  032 

.00050129 

.00050121 

.00000008 

.01551002 

.033 

.0005326^3 

.00053267 

.00000005 

.00870592 

.  034 

.00056507 

.00056507 

.00000000 

.00004424 

.035 

.00059831 

.00059840 

.00000009 

.01552380 

.n36 

.00063265 

.00063266 

.00000001 

.00228879 

.037 

.00066788 

.00066786 

.00000003 

.00403365 

.038 

.00070401 

.00070398 

.00000003 

.00444598 

.039 

.00074102 

.00074102 

.00000000 

.00047502 

.040 

.00077892 

.00077899 

.00000007 

.00946182 

.041 

.00081791 

.00081789 

.00000002 

.00216773 

.o42 

.00085767 

.00085770 

.00000004 

.00412981 

.043 

.00089841 

.00089843 

.00000003 

.00314335 

.044 

.00094002 

.00094008 

.00000006 

.00637750 

.045 

.00098261 

.00098265 

.00000004 

.00382043 

.n46 

.00102607 

.00102613 

.00000006 

.00596453 

.  147 

.00107049 

.00107052 

.00000002 

.00215789 

.  048 

.00111578 

.00111581 

.00000004 

.00335193 
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Table  3.      (cont.) 


TIME 

EXACT 

APPRCX. 

ERROR 

REL.  ERROR 

.049 

.00116202 

.00116202 

.00000000 

.00024096 

.050 

.00120902 

.0012C913 

.00000011 

.00880876 

.051 

.00125708 

.00125715 

.00000006 

.00489227 

.n52 

.00130600 

.00130606 

.00000007 

.00501533 

.053 

.00135586 

.00135588 

.00000002 

.00165947 

.054 

.00140657 

.00140659 

.00000003 

.00204043 

.055 

.00145822 

.00145821 

.00000000 

.00065834 

.056 

.00151060 

.00151071 

.00000011 

.00718918 

.057 

.00156414 

.00156411 

.00000003 

.00163029 

.058 

.00161830 

.00161840 

.00000010 

.00638325 

.059 

.00167350 

.00167358 

.00000008 

.00490589 

.060 

.00172963 

.00172965 

.00000002 

.00116210 

.061 

.00178658 

.00178660 

.00000002 

.00086198 

.062 

.00184436 

.00184443 

.00000007 

.00402308 

.063 

.00190306 

.00190315 

.00000009 

.00471346 

.064 

.00196267 

.00196275 

.00000007 

.00361751 

.065 

.00202321 

.00202322 

.00000000 

.00027184 

.066 

.00208456 

.00208457 

.00000000 

.00039817 

•  067 

.00214672 

.00214680 

.00000008 

.00366140 

.068 

.00220989 

.00220989 

.00000000 

.00030318 

.069 

.00227375 

.00227386 

.00000011 

.00486421 

.070 

.00233863 

.00233870 

.00000007 

.00310866 

.071 

.00240430 

.00240441 

.00000010 

.00424655 

.  072 

.00247087 

.00247098 

.00000010 

.00420904 

.073 

.00253834 

.00253841 

.00000007 

.00293499 

.074 

.00260659 

.00260671 

.00000011 

.00436969 

.075 

.00267574 

.00267586 

.00000013 

.00483605 

.o76 

.00274576 

.00274588 

.00000012 

.00423562 

.077 

.00281658 

.00281675 

.00000017 

.00620612 

.078 

.00288837 

.00288848 

.00000011 

.00366989 

.079 

.00296095 

.00296106 

.00000011 

.00375894 

.080 

.00303429 

.00303449 

.00000020 

.00655508 

.081 

.00310861 

.00310877 

.00000016 

.00521134 

.082 

.00318369 

.00318389 

.00000020 

.00636996 

.083 

.00325976 

.00325987 

.00000011 

.00348492 

.084 

.00333647 

.00333669 

.00000022 

.00658181 

.085 

.003+1415 

.00341435 

.00000020 

.00573788 

.086 

.00349269 

.00349285 

.00000015 

.00440634 

.087 

.00357199 

.00357219 

.00000020 

.00562432 

.088 

.00365214 

.00365237 

.00000022 

.00613339 

.089 

.00373314 

.00373338 

.00000024 

.00636194 

.090 

.00381499 

.00381522 

.00000024 

.00620448 

.091 

.00389758 

.00389790 

.00000032 

.00825641 

.092 

.00398122 

.00398141 

.00000019 

.00467444 

.093 

.00406550 

.00406575 

.00000024 

.00602140 

.094 

.00415062 

.00415091 

.00000029 

.00706401 

.095 

.00423667 

.00423690 

.00000023 

.00546420 

.096 

.00432345 

.00432371 

.00000026 

.00603916 

.097 

.00441106 

.00441134 

.00000028 

.00639755 

.o98 

.00449950 

.00449980 

.00000030 

.00655628 

.099 

.00458876 

.00458907 

.00000030 

.00660091 

ERROR    LC-2     IN       STATEMENT    00C1    +    00    LINES 
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Table  k. 

SOLUTION  OF  D( 2 )Y/DT( 2 ) +2DY/DT+Y= 1 , MODIFIED  SIMPSONS  CONV. 


TIME 

EXACT 

APPROX. 

ERROR 

REL.  ERROR 

.030 

.00044113 

.00044110 

.00000003 

.00672813 

.n4C 

.00077892 

.00077899 

.00000007 

.00884686 

.050 

.00120902 

.00120911 

.00000009 

.00733650 

.060 

.00172963 

.00172961 

.00000002 

.00117366 

.  ")70 

.00233863 

.00233863 

.00000000 

.00004704 

.080 

.00303429 

.00303436 

.00000008 

.00247175 

.090 

.00381499 

.00381503 

.00000004 

.00107733 

.i  on 

.00467886 

.00467886 

.00000000 

.00013037 

•  no 

.00562414 

.00562414 

.00000000 

.00007468 

.1  20 

.00664914 

.00664914 

.00000000 

.00000902 

.130 

.00775208 

.00775220 

.00000012 

.00150411 

.140 

.00893154 

.00893165 

.00000011 

.00122599 

.150 

.01018580 

.01018587 

.OOOOuOO  / 

.00067741 

.160 

.01151319 

.01151325 

.00000006 

.00054720 

.170 

.01291218 

.01291222 

.00000004 

.00030979 

.180 

.01438115 

.01438121 

.00000006 

.00044503 

.190 

.015^1855 

.01591870 

.00000015 

.00096114 

.200 

.01732314 

.01752318 

.00000004 

.00021686 

.710 

.01919310 

.01919315 

.00000005 

.00028135 

.?20 

.02092706 

.02092717 

.00000011 

.00052564 

.730 

.02272367 

.02272379 

.00000012 

.00050608 

.?40 

.02458151 

.02458158 

.00000007 

.00028883 

.  2  5  0 

.02649910 

.02649916 

.00000006 

.00023774 

.260 

.02847497 

.02847516 

.00000019 

.00065672 

.270 

.03050803 

.03050821 

.00000018 

.00059001 

•  280 

.03259673 

.03259699 

.00000026 

.00080069 

.290 

.03473995 

,03474019 

.00000024 

.00069660 

.-^00 

.03693632 

..3693652 

.00000020 

.00054689 

.310 

.03918442 

.03918471 

.00000029 

.00074519 

.320 

.04148330 

.04148351 

.00000021 

.00051105 

.■330 

.04383135 

.04383169 

.00000034 

.00078026 

.^40 

.04622777 

.04622804 

.00000027 

.C0058839 

.350 

.04867106 

.04867137 

.00000031 

.00063898 

.  ->60 

.05116020 

.05116051 

.00000031 

.00059812 

.370 

.05369398 

.05369430 

.00000032 

.00059225 

.380 

.05627125 

.05627162 

.00000037 

.00065398 

.390 

.05889090 

.05889134 

.00000044 

.00074375 

.400 

.06155196 

.06155236 

.00000040 

.00064823 

.410 

.06425318 

.06425360 

.00000042 

.00065522 

.420 

.06699353 

.06699401 

.00000048 

.00071201 

.430 

.06977197 

.06977253 

.00000056 

.00079545 

.440 

.07258757 

.07258813 

.00000056 

.00076735 

.450 

.07543921 

.07543980 

.00000059 

.00078606 

.  V60 

.07832591 

.07832655 

.00000064 

.00082093 

.470 

.08124662 

.08124740 

.00000078 

.00095512 

.480 

.08420055 

.08420137 

.00000082 

.00097862 

.490 

.08718664 

.08718753 

.00000089 

.00102539 

.^00 

.09020405 

.09020495 

.00000090 

.00099330 
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Table  k.     (cont.) 

TIME       EXACT     APPRCX.      ERROR      REL.  ERROR 


.510 

.09325164 

.09325269 

.00000105 

.00112599 

.520 

.09632880 

.09632987 

.00000107 

.00110559 

.530 

.09943434 

.09943558 

.00000124 

.00125007 

.540 

.10256773 

.10256898 

.00000125 

.00121871 

.550 

.10572780 

.10572920 

.00000140 

.00132416 

.  560 

.10891380 

.10891539 

.00000159 

.00145987 

.570 

.11212  509 

.11212671 

.00000162 

.00144481 

.580 

.11536052 

.11536235 

.00000183 

.00158633 

.590 

.11861970 

.11862150 

.00000180 

.00151745 

.600 

.12190140 

.12190338 

.00000198 

.00162426 

.610 

.12520505 

.12520721 

.00000216 

.00172517 

.620 

.12853004 

.12853222 

.00000218 

.00169610 

.630 

.13187535 

.13187768 

.00000233 

.00176682 

.640 

.13524043 

.13524285 

.00000242 

.00178941 

.650 

.13862454 

.13862700 

.00000246 

.00177458 

.660 

.14202682 

.14202943 

.00000261 

.00183768 

.670 

.14544675 

.14544944 

.00000269 

.00184947 

.680 

. 14888344 

.14888634 

.00000290 

.00194783 

.690 

.15233650 

.15233946 

.00000296 

.00194307 

.700 

.15580497 

.15580813 

.00000316 

.00202618 

.710 

.15928841 

.15929170 

.00000329 

.00206544 

.  '20 

.16278616 

.16278954 

.00000338 

.00207634 

.7  30 

.16629754 

.16630101 

.00000347 

.00208662 

.740 

.16982178 

.16982549 

.00000371 

.00218464 

.750 

.17335858 

.17336238 

.00000380 

.00219199 

.760 

.17690710 

.17691108 

.00000398 

.00224977 

.770 

.18046693 

.18047101 

.00000408 

.00226080 

.780 

.18403730 

.18404159 

.00000429 

.00233105 

.79  0 

.18761781 

.18762226 

.00000445 

.00237184 

.800 

.19120791 

.19121247 

.00000456 

.00238484 

.810 

.19480696 

.19481168 

.00000472 

.00242291 

.820 

.19841444 

.19841934 

.00000490 

.00246958 

.830 

.20202989 

.20203492 

.00000503 

.00248973 

.940 

.205  >5265 

.20565792 

.00000527 

.00256257 

.850 

.20928238 

.20928783 

.00000545 

.00260414 

.«6C 

.21291849 

.21292415 

.00000566 

.00265829 

.870 

.21656063 

.21656640 

.00000577 

.00266438 

.R80 

.22020812 

.22021410 

.00000598 

.00271561 

.890 

.22386066 

.22386678 

.00000612 

.00273384 

.900 

.22751769 

.22752398 

.00000629 

.00276462 

.910 

.23117873 

.23118525 

.00000652 

.00282033 

.920 

.23484346 

.23485015 

.00000669 

.00284871 

.930 

.23851134 

.23851824 

.00000690 

.00289294 

.940 

.24218202 

.24218910 

.00000708 

.00292342 

.950 

.24585500 

.24586231 

.00000731 

.00297330 

.960 

.24952992 

.24953746 

.00000754 

.00302168 

.970 

.25320644 

.25321415 

.00000771 

.00304495 

.980 

.25688402 

.25689198 

.00000796 

.00309867 

.Q90 

.26056237 

.26057057 

.00000820 

.C33147C4 
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Table  5. 

SOLUTION  OF  D(2)Y/DT(2)+2DY/DT+Y=1» TRAPEZOIDAL  CONVOLUTION 


TIME 

EXACT 

APPROX. 

ERROR 

REL.  ERROR 

.0200 

.00019733 

.00019703 

.00000029 

.14767441 

.0300 

.00044113 

.00044063 

.00000050 

.11328330 

•  0400 

.00077892 

.00077837 

.00000055 

.07003415 

.n500 

.00120902 

.00120836 

.00000066 

.05487071 

.0600 

.00172963 

.00172872 

.00000090 

.05223673 

.0700 

.00233863 

.00233762 

.00000101 

.04314923 

.0800 

.00303429 

.00303324 

.00000105 

.03466052 

.0900 

.00381499 

.00381379 

.00000120 

.03141033 

.1000 

.00467886 

.00467752 

.00000134 

.02864161 

.]  100 

.00562414 

.00562269 

.00000145 

.02581195 

.1200 

.00664914 

.00664760 

.00000154 

.02320601 

.1300 

.00775208 

.00775056 

.00000152 

.01957668 

.1400 

.00893154 

.00892993 

.00000161 

.01802937 

.15  0  0 

.01018580 

.01018407 

.00000173 

.01699425 

.1600 

.01151319 

.01151138 

.00000181 

.01575584 

.1700 

.01291218 

.01291027 

.00000191 

.01479224 

.18  0  0 

.01438115 

.01437920 

.00000196 

.01359418 

.1  900 

.01591855 

.01591662 

.00000193 

.01213050 

.?000 

.01752314 

.01752103 

.00000211 

.01202410 

.2100 

.01919310 

.01919095 

.00000215 

.01119673 

.2200 

.02092706 

.02092491 

.00000215 

.01026*00 

.2300 

.02272367 

.02272147 

.00000220 

.00966393 

.  >400 

.02458151 

.02457922 

.00000229 

.00930781 

.2500 

.02649910 

.02649676 

.00000234 

.00883426 

.?600 

.02847497 

.02847271 

.00000226 

.00793328 

.7700 

.03050803 

.03050573 

.00000230 

.00754883 

.2800 

.03259673 

.03259447 

.00000226 

.00692094 

.?900 

.03473995 

.03473764 

.00000231 

.00664077 

.3000 

.03693632 

.03693394 

.00000238 

.00643269 

.3100 

.03918442 

.03918211 

.00000231 

.00590286 

.3200 

.04148330 

.04148088 

.00000242 

.00582885 

.3300 

.04383135 

.04382904 

.00000231 

.00527020 

.3400 

.04622777 

.04622537 

.00000240 

.00519168 

.3500 

.04867106 

.04866868 

.00000238 

.00488*97 

.3600 

.05:16020 

.05115780 

.00000240 

.00469310 

.3700 

.05369398 

.05369157 

.00000241 

.00448095 

.3800 

.05627125 

.05626887 

.00U00238 

.00422774 

.^900 

.05889090 

.05888858 

.00000232 

.00393270 

.4000 

.06155196 

.06154960 

.00000236 

.00382766 

.4100 

.06425318 

.06425085 

.00000233 

.00362472 

.4200 

.06699353 

.06699127 

.00000226 

.00337794 

.4300 

.06977197 

.06976979 

.00000218 

.00312160 

.4400 

.07258757 

.07258541 

.00000216 

.00298123 

.4500 

.07543921 

.07543709 

.00000212 

.00281154 

.4600 

.07832591 

.07832384 

.00000207 

.00264025 

.4700 

.08124662 

.08124469 

.00000193 

.00237179 

.4800 

.08420055 

.08419867 

.00000188 

.00222920 

.4900 

.08718664 

.08718484 

.00000180 

.00206224 
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Table  5.  (cont.) 


TIME 

EXACT 

APPROX. 

ERROR 

REL.  ERROR 

.  iOOO 

.09020405 

.09020227 

.00000178 

.00197663 

.5100 

.09325164 

.09325002 

.00000162 

.00173723 

.5200 

.09632880 

.09632720 

.00000160 

.00166305 

.6300 

.09943434 

.09943291 

.00000143 

.00143713 

.5400 

.10256773 

.10256630 

.00000143 

.00139420 

.5500 

.10572780 

.10572650 

.00000130 

.00122957 

.5600 

.10891380 

.10691266 

.00000114 

.00104670 

.5700 

.11212509 

.11212395 

.U0u00114 

.00101672 

.5800 

.11536052 

.11535956 

.00000096 

.uu063217 

.5900 

.11861970 

.11861868 

.00000102 

.00085989 

.6000 

.12190140 

.12190053 

.00000087 

.00071369 

,6100 

.12520505 

.12520434 

.00000071 

.00056707 

.6200 

.12153004 

.12852935 

.00000069 

.00053684 

.6300 

.13187535 

.13187480 

.00000055 

.00041706 

.6400 

.13524043 

.13523995 

.00000048 

.00035492 

.6500 

.13862454 

.13862409 

.00000045 

.00032462 

.6600 

.14202682 

.14202650 

.00000032 

.00022531 

.6700 

.14544675 

.14544649 

.00000026 

.00017876 

.6800 

.14888344 

.14888337 

.00000007 

.00004702 

.6900 

.15233650 

.15233646 

.00000004 

.0u002626 

.7000 

.15580497 

.15580510 

.00000013 

.00008344 

.  7  1 C  0 

.15928841 

.15928863 

.00000022 

.00013811 

.7200 

.16278616 

.16278642 

.00^00026 

.00015972 

.7300 

.16629754 

.16629784 

.00000030 

.00018040 

.7400 

.16982178 

.16982227 

.00000049 

.00028854 

.7500 

.17335858 

.17335911 

.00000053 

.00030572 

.7600 

.17690710 

.17690775 

.00000065 

.00036742 

.7700 

.18046693 

.18046761 

.00000068 

.30037680 

.7800 

.18403730 

.16403812 

.00000082 

.00044556 

.7900 

.18761781 

.18761872 

.000C0O91 

.00048503 

•  8000 

.19120791 

.19120885 

.00000094 

.00049161 

.8100 

.1940C696 

.19480797 

.00000101 

.00051846 

.8200 

.19841444 

.19841554 

.00000110 

.00055440 

.8300 

.20202989 

.20203104 

.00000115 

.00056922 

.8400 

.20565265 

.20565396 

.00000131 

.00063700 

.850.0 

.20928238 

.20928378 

.00000140 

.00066895 

.8600 

.21291849 

.21292001 

.00000152 

.00071389 

.R7P0 

.21656063 

.21656217 

.00000154 

.00071112 

.8800 

.22020812 

.22020977 

.00000165 

.00074929 

.  }900 

.22386066 

.22386235 

.00000169 

.00075493 

.9000 

.22751769 

.22751945 

.00000176 

.00077357 

.9100 

.23117873 

.23118061 

.00  000188 

-5  132  2 

.9200 

.23484346 

.23484539 

.00000193 

.00082182 

.9300 

.23851134 

.23851337 

.00000203 

.00085111 

.9400 

.24218202 

.24218411 

.00000209 

.00086299 

.9500 

.24585500 

.24565719 

.00000219 

.00089077 

.9600 

.24952992 

.24953220 

.00000228 

.00091372 

.9700 

.25320644 

.25320874 

•00U00230 

.00090835 

.9800 

.25688402 

.25688642 

.00000240 

.00093427 

.9900 

.26056237 

.26056485 

.00000248 

.00095179 

ERROR  LC-2  IN   STATEMENT  0001  +  00  LINES 


36 
Table  6. 
SOLUTION    OF    D(2)Y/DT(2)+2DY/OT+Y=l»TRAPEZCIDAL    CONVOLUTION 


TTME 

EXACT 

APPROX. 

ERROR 

REL.  ERROR 

,o020 

.00000190 

.00000200 

.00000010 

5.33830050 

.0030 

.00  100439 

.00000449 

.00000010 

2.37822030 

.0040 

.00000787 

.00000798 

.00000011 

1.40370630 

.0050 

.00001244 

.00001246 

.00000002 

.16150352 

.0060 

.00001789 

.00001793 

.00000004 

.19659624 

.0070 

.00002433 

.00002438 

.00000006 

.23148982 

.0080 

.00003174 

.00003183 

.00000008 

.26630229 

.0090 

.00004024 

.00004026 

.00000002 

.04957949 

.0100 

.00004960 

.00004967 

.00000007 

.13162235 

.0110 

.00006004 

.00006006 

.00000002 

.03227010 

.0120 

.00007134 

.00007142 

.00000009 

.12145701 

.0130 

.00008371 

.00008377 

.00000006 

.07662744 

.0140 

.00009703 

.00009709 

.00000005 

.05553930 

.0150 

.00011132 

.00011138 

.00000006 

.05294646 

.0160 

.00012656 

.00012664 

.00000008 

.06189062 

.0170 

.00014286 

.00014287 

.00000001 

.00917001 

.0180 

.00016000 

.00016007 

.00000007 

.04143750 

.0190 

.00017819 

.00017823 

.00000004 

.02019182 

.0200 

.00019733 

.00019735 

.00000002 

.01226397 

.0210 

.00021740 

.00021743 

.00000003 

.0151885^ 

.0220 

.00023841 

.00023848 

.00000006 

.02663026 

.0230 

.00026046 

.00026048 

.00000001 

.00524837 

.0240 

.00028344 

.00028343 

.00000001 

.00436068 

.0250 

.00030735 

.00030734 

.00000002 

.00505284 

.0260 

.00033219 

.00033219 

.00000000 

.00247753 

.0270 

.00035805 

.00035800 

.00000005 

.01284748 

.0280 

.00038472 

.00038475 

.00000003 

.00827868 

.0290 

.00041242 

.00041245 

.00000003 

.00821250 

.  1300 

.00044113 

.00044110 

.00000004 

.00835802 

.0310 

.00047065 

.00047068 

.00000003 

.00546474 

.032  0 

.00050129 

.00050120 

.00000009 

.01720764 

.0330 

.00053263 

.00053266 

.00000004 

.00701055 

.n340 

.00056507 

.00056506 

.00000000 

.00167413 

.0350 

.00059831 

.00059839 

.00000008 

.01379058 

.0360 

.00063265 

.00063265 

.00000000 

.00052952 

.^370 

.00066788 

.00066784 

.00000004 

.00581540 

.0380 

.00070401 

.00070396 

.00000004 

.00624426 

.0390 

.00074102 

.00074101 

.00000000 

.00133735 

.n4C0 

.00077892 

.00077898 

.00000006 

.00762851 

.0410 

.00081791 

.00081787 

.00000003 

.00400658 

.0420 

.00185767 

.00085769 

.00000002 

.00227478 

.0^30 

.00089841 

.00089842 

.00000001 

.00127782 

.0440 

.00094002 

.00094007 

.00000004 

.00449244 

.0450 

.00098261 

.00098263 

.00000002 

.00192039 

.0460 

.00102607 

.00102611 

.00000004 

.00405432 

.0470 

.00107049 

.00107049 

.00000000 

.00023354 

.0480 

.00111578 

.00111579 

.00000002 

.00139813 
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Table  6.     (cont.) 


TIME 

EXACT 

APPRCX. 

ERROR 

REL.  ERROR 

•  0490 

.00116202 

.00116200 

.00000003 

.00222887 

.0500 

.00120902 

.00120911 

.00000008 

.00679060 

.0510 

.00125708 

.00125712 

.00000004 

.00284786 

.0520 

.00130600 

.00130604 

.00000004 

.00293263 

.0530 

.00135586 

.00135585 

.00000000 

.00046465 

•  0540 

.00140657 

.00140656 

.00000000 

.00013508 

.0550 

.00145822 

.00145817 

.00000004 

.00289395 

.0560 

.001.51060 

.00151068 

.00000007 

.00489208 

.0570 

.00156414 

.00156407 

.00000006 

.00398942 

.0580 

.00161830 

.00161836 

.00000006 

.00395477 

.0590 

.00167350 

.00167354 

.00000004 

.00241411 

.0600 

.00172963 

.00172960 

.00000002 

.00139336 

.0610 

.00178658 

.00178655 

.00000003 

.00175755 

.0620 

.001&4436 

.00184438 

.00000002 

.00133922 

.0630 

.00190306 

.00190310 

.00000004 

.00197051 

.0640 

.00196267 

.00196269 

.00000002 

.00082031 

.0650 

.00202321 

.00202316 

.00000005 

.00258994 

.0660 

.00208456 

.00208451 

.00000005 

.00252811 

.0670 

.00214672 

.00214673 

.00000001 

.00067079 

.0680 

.00220989 

.00220983 

.00000006 

.00275580 

.0690 

.00227375 

.00227379 

.00000004 

.00173722 

.0700 

.00233863 

.00233863 

.00000000 

.00009407 

.0710 

.00240430 

.00240433 

.00000002 

.00096910 

.0720 

.00247087 

.00247089 

.00000002 

.00085395 

.0730 

.00253834 

.00253832 

.00000001 

.00049639 

.0740 

.00260659 

.00260662 

.00000002 

.00086320 

.0750 

.00267574 

.00267577 

.00000003 

.00125199 

.0760 

.00274576 

.00274578 

.00000002 

.00056815 

.077C 

.00281658 

.00281665 

.00000007 

.00245688 

.0780 

.00288837 

.00288837 

.00000000 

.00015926 

.0790 

.00296095 

.00296094 

.00000000 

.00015536 

.o800 

.00303429 

.00303437 

.00000008 

.00255085 

.0810 

.00310861 

.00310864 

.00000003 

.00111626 

.0820 

.00318369 

.00318376 

.00000007 

.00218300 

.  1830 

.00325976 

.00325973 

.00000003 

.00079761 

.0840 

.00333647 

.00333654 

.00000007 

.00220293 

.0850 

.00341415 

.00341419 

.00000004 

.00126532 

.0860 

.00349269 

.00349269 

.00000000 

.00016320 

.0870 

.00357199 

.00357202 

.00000003 

.00095465 

.0880 

.00365214 

.00365219 

.00000005 

.00136358 

.0890 

.00373314 

.00373320 

.00000006 

.00148936 

.0900 

.00381499 

.00381503 

.00000005 

.00122936 

.0910 

.00389758 

.00389770 

.00000012 

.00317890 

.0920 

.00398122 

.00398120 

.00000002 

.00050738 

.0930 

.00406550 

.00406553 

.00000003 

.00073546 

,o940 

.00415062 

.00415069 

.00000007 

.00167204 

.0950 

.0C'f23667 

.00423667 

.00000000 

.00003541 

.0960 

. 00*32345 

.00432347 

.00000002 

.00043021 

.0970 

.00441106 

.00441109 

.00000003 

.00068011 

.0980 

.00449950 

.00449953 

.00000003 

.00073341 

.0990 

.00458876 

.00458879 

.00000003 

.00067120 

ERROR  LC-2  IN   STATEMENT  0001  +  Oo 

LINES 

38 

Table  7. 

SOLUTION  OF  DO)  Y/ DT  (  3  ) +  3D  (  2  )  Y/DT  (  2  )+3DY/DT+Y=l  ,MOD.  SIMP.  CCNV. 
TIMl £XmCT     aPPRCX. lKKCR klL.  iRROR 


.0400 

.00001029 

.00001068 

.00000039 

3.79652210 

.0500 

.00001999 

.00002049 

.00000050 

2.50555110 

.0600 

.00003445 

.00003494 

.00000049 

1.41191240 

.0700 

.00005426 

.00005487 

.00000060 

1.11429520 

.0800 

.00008032 

.00008110 

.00000078 

.97599218 

.0900 

.00011357 

.00011442 

.00000086 

.75470170 

.1000 

.00015467 

.00015560 

.00000093 

.60307183 

.1100 

.00020434 

.00020538 

.00000104 

.50750894 

.1100 

.00026331 

.00026447 

.00000116 

.44023356 

.1300 

.00031217 

.00033357 

.00000140 

.42027799 

.1400 

.0004x183 

.00041334 

.00000151 

.36725716 

.1500 

.00050284 

.00050443 

.00000160 

.31791741 

.1600 

.00060575 

.00060747 

.00000172 

.28409457 

.1700 

.00072122 

.00072305 

.00000183 

.25347500 

.1800 

.00084977 

.00085176 

.00000198 

.23358148 

.1900 

.00099194 

.00099415 

.00000221 

.22277445 

.2000 

.00114852 

.00115076 

.00000223 

.19429285 

•  2100 

.00131972 

.00132210 

.00000239 

.16091771 

.2200 

.00150611 

.00150869 

.00000258 

.17150199 

.2300 

.0017U826 

.00171099 

.00000273 

.15968236 

.2400 

.00192663 

.00192946 

.00000283 

.14707569 

.2500 

.00216158 

.00216455 

.00000297 

.13756949 

.2600 

.00241343 

.00241668 

.00000325 

.13468825 

.2700 

.00268285 

.00268625 

.00000340 

.12669754 

.2800 

.00297001 

.00297365 

.00000364 

.12256682 

.2900 

.00327547 

.00327925 

.00000376 

.11546147 

.3000 

.00359950 

.00360340 

.00000391 

.10854294 

.^100 

.00394229 

.00394645 

.00000416 

.10562133 

.3200 

.00430447 

.00430872 

.00000425 

.09880199 

.3300 

.00468595 

.00469051 

.00000456 

.09731220 

.•3400 

.00508744 

.00509211 

.00000467 

.09174353 

.^500 

.00550891 

.00551380 

.00000489 

.06873983 

.3600 

.00595077 

•005955o4 

.00000507 

.08525786 

.3700 

.00641^21 

.00641fa4v 

.00000526 

.0o^6300 

.3800 

.00689645 

.006*0196 

.00000551 

.07987584 

.3900 

.00740072 

.00740649 

.00000577 

.07795727 

.40  10 

.00792635 

.00793228 

.00000593 

.07478345 

.4100 

.00847337 

.00847953 

.00000615 

.07262633 

.4200 

.00904200 

.00904841 

.00000641 

.07091795 

.4300 

.00963240 

.00963910 

.00000670 

.06960364 

.4400 

.01024484 

.01025176 

.00  .00692 

.06751689 

.4500 

.01087936 

.01088653 

.00000717 

.06593221 

.4600 

.01153610 

.01154354 

.00000745 

.06453656 

.4700 

.01221512 

.01222292 

.00000781 

.06391262 

.4800 

.01291670 

.01292479 

.00000809 

.06260113 

.4900 

.01364084 

.01364923 

.00000839 

.06153582 

.5000 

.01438771 

.01439635 

.00000864 

.06001648 
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Table  7.  (cont.) 

TIME EXACT     APPROX.      ERROR      REL.  ERROR 


.5100 

.01515719 

.01516622 

.00000904 

.05962189 

.^200 

.01594962 

.01595892 

.00000931 

.05833995 

.^300 

.01676477 

.01677451 

.00000975 

.05814576 

.  5400 

.0176^304 

.01761305 

.00001001 

.05684245 

.  5500 

.0184o414 

.01847456 

.00001042 

.05645538 

.5600 

.01934822 

.01935910 

.00001068 

.05624290 

.5700 

.02025548 

.02026669 

.00001121 

.05532330 

.  5800 

.02118561 

.02119733 

.00001172 

.05533^43 

.5900 

.02213904 

.02215106 

.00001202 

.05427968 

.6000 

.02311530 

.02312786 

.00001256 

.05432333 

.6100 

.02411462 

.02412773 

.00001311 

.05435707 

.6200 

.02513712 

.02515066 

.00001354 

.05385661 

.6300 

.02618250 

.02619663 

.00001413 

.05395207 

.6400 

.02725094 

.02726560 

.00001466 

.05381099 

.6500 

.02834237 

.02835756 

.00001519 

.05358762 

.6600 

.02945660 

.02947245 

.00001585 

.05379440 

.6700 

.03059376 

.03061022 

.00001646 

.05380182 

.6800 

.03175359 

.03177083 

.00001724 

.05428048 

.6900 

.03293631 

.03295420 

.00001789 

.05432910 

.7000 

.03414156 

.03416029 

.00001873 

.w5485104 

.7100 

.03536949 

.03538901 

.00001952 

.055l83i;> 

.7200 

.03661997 

.03664029 

.00002032 

.05549158 

.7300 

.03789289 

.03791406 

.00002117 

.05586272 

.7400 

.03918798 

.03921021 

.00002223 

.05673168 

.7500 

.04050548 

.04052865 

.00002317 

.u5720461 

.7600 

.04184503 

.04186927 

.00002424 

.05792803 

.7700 

.04320670 

.04323197 

.00002527 

.05847473 

.7800 

.04459019 

.04461663 

.00002644 

.05929331 

.7900 

.04599554 

.04602315 

.00002761 

.06002756 

.8000 

.04742263 

.04745142 

.00002679 

.06070098 

.81  ">0 

.048tS7127 

.04890131 

.00003004 

.06146351 

.82U0 

.05034132 

.05037269 

.00003137 

.06232058 

.8300 

.05183271 

.051d6545 

.00003274 

.06316089 

.8400 

.05334517 

.05337945 

.00003426 

.06425137 

.8500 

.05487873 

.05491455 

.00003562 

.06526572 

.8600 

.05643315 

.05647062 

.00003747 

.06639006 

.R700 

.05800840 

.05804751 

.00003911 

.06741610 

.R800 

.05960417 

.05964507 

.00004090 

.06861936 

.8900 

.06122044 

.06126315 

.00004271 

.06976592 

.9000 

.06285697 

.06290160 

.00004463 

.07100088 

.9100 

.06451356 

.06456026 

.00004670 

.07239098 

.9100 

.06619019 

.06623698 

.00004879 

.07371032 

.9300 

.067tS°>658 

.06793759 

.00005101 

.07513414 

.9400 

.06960264 

.06965592 

.00005328 

.07654307 

.9500 

.07133810 

.07139380 

.00005570 

.07808310 

.9600 

.07309287 

.07315108 

.00005821 

.07963841 

.9700 

.07486682 

.07492758 

.00006076 

.08115210 

.Q800 

.07665963 

.07672312 

.00006349 

.08282065 

.9900 

.07847121 

.07853753 

.00006632 

.08451635 
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Table  8. 
."LUTICN    OF    D(3)Y/DT(3)+3D(2 )Y/DT(2 )+3DY/DT+Y=l ,MOD.    SIMP.    CCNV. 


TTMF 

EXACT 

APPROX. 

ERROR 

REL.  ERROR 

.n040 

-.00000010 

.00000001 

.00000011- 

110.91657000 

.0050 

-.00000000 

.00000002 

.00000002 

.0060 

-.00000000 

.00000004 

.00000004 

.0070 

-.oooooooc 

.00000006 

.00000006 

.0080 

-.00000000 

.00000009 

.00000009 

.0090 

.00000010 

.00000012 

.00000002 

22.34238000 

.0100 

.00000010 

.00000017 

.00000007 

68.79092400 

.0110 

.00000020 

.00000022 

.00000002 

11.10989600 

.0120 

.0000020 

.00000029 

.00000009 

45.52710000 

.0130 

.00000030 

.00000036 

.00000007 

22.71997600 

.0140 

.00000040 

.00000045 

.00000006 

14.45889700 

.0150 

.00000049 

.00000056 

.00000006 

12.64812900 

.0160 

.00000059 

.00000068 

.00000008 

13.84765100 

.0170 

.00000079 

.00000081 

.00000002 

2.16957980 

.0180 

.00000089 

.00000096 

.00000007 

7.93652320 

.ol90 

.00000109 

.00000113 

.00000004 

3.72433190 

.0200 

.00000129 

.00000132 

.00000003 

2.25215740 

.0210 

.00000148 

.00000152 

.00000004 

2.56267150 

.0220 

.00000168 

.00000175 

.00000007 

4.09320720 

.0230 

.00000198 

.00000200 

.00000002 

.96763256 

.0240 

.00000227 

.00000227 

.00000000 

.29443459 

.0250 

.00000257 

.00000256 

.00000000 

.37543668 

.0260 

.00000286 

.00000288 

.00000001 

.50181423 

.0270 

.00000326 

.00000322 

.00000004 

1.21358990 

.0280 

.00000355 

.00000359 

.00000004 

1.08279460 

.0290 

.00000394 

.00000398 

.00000004 

1.03136060 

.0300 

.00000443 

.00000440 

.00000003 

.67516227 

.0310 

.00000482 

.00000485 

.00000003 

.68225284 

.0320 

.00000541 

.00000533 

.00000008 

1.45832650 

.0330 

.00000580 

.00000585 

.00000004 

.77213281 

.0340 

.00000639 

.00000639 

.00000000 

.02987950 

.0350 

.00000687 

.00000696 

.00000009 

1.31191130 

.0360 

.00000756 

.00000757 

.00000001 

.14505638 

.0370 

.00000825 

.00000821 

.00000003 

.37992258 

.0380 

.00000893 

.00000889 

.00000004 

.40966417 

.03  >0 

.00000961 

.00000960 

.00000000 

.02776367 

.0400 

.00001029 

.00001035 

.00000007 

.64553346 

.0410 

.00001117 

.00001114 

.00000003 

.23463081 

.042  0 

.00001194 

.00001197 

.00000003 

.21556102 

.0430 

.00001282 

.00001283 

.00^00002 

.13443251 

.0440 

.00001369 

.00001374 

.00000005 

.34595531 

•  0450 

.00001466 

.00001469 

.00000002 

.15942999 

.046  0 

.00001563 

.00001568 

.00000005 

.29029782 

.n470 

.00001670 

.00001671 

.00000000 

.03160006 

.0480 

.00001777 

.00001778 

.00000002 

.09930878 

.0490 

.00001893 

.00001890 

.00000002 

.13085031 

•  0  "5  0  0 

.00001999 

.00002007 

.00000008 

.41228092 
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Table  8.      (cont.) 
ME EXACT  APPRCX.  ERROR  REL.    ERROR 


.0510 

.00002125 

.00002128 

.00000004 

.16526487 

.0520 

.00002250 

.00002254 

.00000004 

.16270668 

.0530 

.00002386 

.00002385 

.00000000 

.03746039 

.0540 

.00002521 

.00002520 

.00000000 

.02181277 

.0550 

.00002666 

.00002661 

.00000005 

.17528097 

.0560 

.00002800 

.000028C7 

.00000007 

.24474883 

.0570 

.00002964 

.00002957 

.00000007 

.23244580 

.0580 

•00003108 

.00003114 

.00000006 

.18217037 

.0590 

.00003272 

.00003275 

.00000003 

.09696853 

.0600 

.00003445 

.00003442 

.00000003 

.09883657 

.0610 

.00003618 

.00003614 

.00000004 

.11764966 

.0620 

.00003790 

.00003792 

.00000001 

.03244070 

.0630 

.00003973 

.00003975 

.00000002 

.05954347 

.^640 

.00004164 

.00004164 

.oooooooo 

.00200290 

.06  JO 

.00004366 

.00004359 

.00000007 

.15862745 

.0660 

.00004567 

.00004560 

.00000007 

.15590040 

.0670 

.00004767 

.00004767 

.oooooooo 

.01213101 

.0680 

.00004988 

.00004979 

.00000008 

.16679899 

•  0690 

.000^5197 

.00005198 

.00000002 

.02977432 

.0700 

.00005426 

.00005423 

.00000003 

•05229162 

.0710 

.00005655 

.00005655 

.OOOOOOOO 

.00927074 

.0720 

.OO0U5694 

.00005893 

.oooooooo 

.01634463 

.0730 

.00006142 

.00006137 

.00000005 

.07448608 

.0740 

.00006389 

.00006388 

.00000001 

.02080578 

.0750 

.00006646 

.00006645 

•oooooooo 

.00720473 

.0760 

.00006912 

.00006909 

.00000003 

.03666939 

.0770 

.0000"'178 

.00007180 

.00000003 

.03526422 

.0780 

.0000/463 

.00007458 

.00000005 

.06932149 

.n790 

.00007748 

.00007742 

.00000006 

.07112563 

.0800 

.000^8032 

.00008034 

.00000002 

.02947376 

.0810 

.00008335 

.00008333 

.00000002 

.027021513 

.0820 

•00008638 

.00008638 

.oooooooo 

.  .981645 

.0830 

.00008961 

.00008951 

.00000009 

.10185089 

.0840 

.00009271 

.00009272 

.oooooooo 

.00475445 

.o850 

.00009603 

.00009599 

.00000003 

.03169585 

.0860 

.00009943 

.00009935 

.00000008 

.08442388 

.0870 

.00010282 

.00010277 

.00000005 

.04751888 

.0880 

.00010631 

.00010628 

.00000004 

.03605331 

.0890 

.00010989 

.00010986 

.00000004 

.03421482 

.0900 

.00011357 

.00011351 

.00000005 

.04554204 

.0910 

01172  3 

.00011725 

•00000002 

.01667665 

.0920 

.00012119 

.00012106 

.00000013 

.10796659 

.093  0 

.00012505 

.00012496 

.00000009 

.06964659 

.0940 

.00012899 

.00012893 

.00000005 

.04197303 

.0950 

.00013312 

.00013299 

.00000013 

.09898484 

.0960 

.00013725 

.00013713 

.00000012 

.08659595 

.0970 

.00014146 

.00014135 

.00000011 

.08126478 

.0980 

.00014577 

.00014565 

.00000012 

.08223768 

.0990 

.00015017 

.00015004 

.00000013 

.08688733 

ERROR  LC-2  IN   STATEMENT  OOol  +  00 

LINES 

Table  9. 
SOLUTION    OF    D(3)Y/DT(3)+3D(2)Y/DT(2 ) +3DY/DT+Y=1 » TRAP.    CCNV. 


TIME 

EXACT 

APPRCX. 

ERROR 

REL.  ERROR 

.0300 

.00000443 

.00000463 

.00000020 

4.42164560 

.0400 

.00001029 

.00001065 

.00000036 

3.48976660 

.0500 

.00001999 

.00002043 

.00000044 

2.19696700 

.0600 

.00003445 

.00003484 

.00000039 

1.12106020 

.0700 

.00005426 

.00005472 

.00000046 

.84381070 

.0800 

.00008032 

.00008090 

.00000058 

.72564075 

.0900 

.00011357 

.00011416 

.00000059 

.52331080 

.1000 

.00015467 

.00015527 

.00000060 

.38869256 

.1100 

.00020434 

.00020497 

.00000063 

. J0825102 

.1200 

.00026331 

.00026398 

.00000067 

.25442393 

.1300 

.00033217 

.00033299 

.00000082 

.24640172 

.1400 

.00041183 

.00041267 

.00000084 

.2u4i2618 

.1500 

.00050284 

.00050366 

.00000083 

.16441974 

.1600 

.00060575 

.00060659 

.00000084 

.13925240 

.1  700 

.00072122 

.00072206 

.00000084 

.11646467 

.1  800 

.00084977 

.00085065 

.00000088 

.10366781 

.19  0  0 

.00099194 

.00099292 

.00000099 

.09930883 

.7000 

.00114852 

.00114941 

.00000088 

.07673327 

.2100 

.00131972 

.00132062 

.00000091 

.06875722 

.2200 

.00150611 

.00150707 

.00000097 

.06429167 

•  2300 

.0017C826 

.00170923 

•0000009b 

.05722210 

.2400 

.00192663 

.00192756 

.00000094 

.04660527 

.2500 

.00216158 

.00216250 

.OOOOOOvj 

.04298717 

.2600 

.00141343 

.00241448 

.00000102 

.U4364751 

.2700 

.00268285 

.00268389 

.00000105 

.038*5117 

.2800 

.00297001 

.00297113 

.00000112 

.03766662 

.2900 

.00327547 

.00327657 

.00000110 

.03358607 

.^000 

.00359950 

.00360055 

.00000105 

.02928743 

•  3100 

.00394229 

.00394343 

.00000114 

.02880051 

.?200 

.00430447 

.00430551 

.00000104 

.02427478 

.3300 

.00468595 

.00468712 

.00000117 

.02492772 

.3400 

.00508744 

.00508853 

.00000109 

.02136633 

.3500 

.00550891 

.00551003 

.00000111 

.02023630 

.3600 

.00595077 

.00595187 

.00000110 

.01650517 

.3700 

.00641321 

.00641431 

.00000110 

.U1714274 

.3800 

.00689645 

.00689758 

.00000112 

.01628083 

.3900 

.00740072 

.00740189 

.00000117 

.01577954 

.^000 

.00792635 

.00792746 

.00000110 

.01391434 

.4100 

.00847337 

.00847447 

.00000110 

.01295942 

.4200 

.00904200 

.00904312 

.00000112 

.01235125 

.4300 

.00963240 

.00963356 

.00000116 

.01204996 

.4400 

.01024484 

.01024596 

.00000112 

.01090305 

.4500 

.01087936 

.01088047 

.00000111 

.01020281 

.4600 

.01153610 

.01153721 

.00000111 

.00963064 

.4700 

.01221512 

.01221631 

.00000119 

.00977477 

.4800 

.01291670 

.01291789 

.00000119 

.00918191 

.4900 

.01364084 

.01364203 

.00000120 

.oo877512 

.5000 

.01438771 

.01438804 

.00U00113 

.00785392 
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Table  9.  (cont.) 
TIME EXACT     APPRCX.      ERROR        L.  ERROR 


.5100 

.01515719 

.01515840 

.00000121 

.00800280 

.5200 

.01594962 

.01595077 

.00000115 

.00721020 

.5300 

.01676477 

.01676602 

.00000125 

.00745015 

.5400 

.01760304 

.01760419 

.00000115 

.00652728 

.5500 

.01846414 

.01846533 

.00000119 

.00646659 

.  c  6  0  C 

.01934822 

.01934948 

.00000126 

.00653290 

.5700 

.02025548 

.02025666 

.00000118 

.00584533 

.  5800 

.02118561 

.02118689 

.00000128 

.00604656 

.5900 

.02213904 

.02214017 

.00000113 

.00512217 

.6000 

.02311530 

.02311651 

.00000121 

.UU524761 

.6100 

.02411462 

.02411590 

.000001^8 

.00531213 

.6200 

.02513712 

.02513832 

.00000120 

.00478575 

.6300 

.02618250 

.02618376 

.00000126 

.00480092 

.6400 

.02725094 

.02725217 

.00000123 

.00452461 

.6500 

.02834237 

.02834353 

.00000116 

,uu410693 

.6600 

.02945660 

.02945780 

.00000120 

.00406700 

.6700 

.03059376 

.03059492 

.00000116 

.00377528 

.6800 

.03175359 

.03175483 

.00000124 

.00390507 

.6900 

.03293631 

.03293748 

.00000117 

.00355535 

.7000 

•  03U4156 

.03414280 

.00000124 

.00363194 

.7100 

.03336949 

.03537072 

.00000123 

.00346909 

.7200 

.03661997 

.03662115 

.00000118 

•00322775 

.7300 

.03789289 

.0378y402 

.00000113 

.0u297681 

.7400 

.03918798 

.03918923 

.00000123 

.uo318465 

.7500 

.04050548 

.04050670 

.00000122 

.00299959 

.7600 

.04184503 

.04184631 

.00000128 

.00306608 

.7700 

.04320670 

.04320798 

.00000128 

.00296250 

.7800 

.04459019 

.04459159 

.00000140 

.00314195 

.7900 

.04599554 

.04599703 

.00000149 

.00324597 

.8000 

.04742263 

.04742419 

.00000156 

.00328324 

•  filOO 

.04887127 

.04887292 

.00000165 

.00338236 

.8200 

.05034132 

.05034312 

.00000180 

.00357758 

.8300 

.05183271 

.05183464 

.00000193 

.00373123 

.8400 

.05334517 

.05334736 

.00000219 

•0041Q346 

.8500 

.05487873 

.05488113 

.00000240 

. J0436781 

.8600 

.05643315 

.05643581 

.00000266 

.u^47o623 

.8700 

.058uU840 

.05801125 

.00000285 

.00491308 

.8800 

.05960417 

.05960731 

.00000314 

.00526809 

.8900 

.06122044 

.06122383 

.00000339 

.00554227 

.9000 

.06285697 

.06286066 

.00000369 

.00587206 

.Q100 

.06451356 

.06451764 

.00000408 

.00631650 

.9200 

.06619019 

.06619459 

.00000440 

.00664298 

.9300 

.06788658 

.06789135 

.00000477 

.00702790 

.9400 

.06960264 

.06960775 

.00000511 

.00734742 

.9500 

.07133810 

.07134362 

.00000552 

.00774341 

.9600 

.07309287 

.07309879 

.00000592 

.00809792 

.  )700 

.07486682 

.07487307 

.00000625 

.00835083 

.9800 

.07665963 

.07666630 

.00000667 

.00869558 

.9900 

.07847121 

.07647827 

.00000706 

.00900075 
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Table  10. 

SOLUTION  OF  D(3)Y/DT(3)+3D(2)Y/DT(2 ) +  3DY/DT  +  Y  =  1 » TRAP.  CCNV. 


TIME 

EXACT 

APPROX. 

ERROR 

REL.  ERROR 

.0030 

-.00000010 

.00000000 

.00000011- 

•104.72249000 

.o040 

-.00000010 

.00000001 

.00000011- 

110.91334000 

.0050 

-.00000000 

.00000002 

.00000002 

.0060 

-.00000000 

.00000004 

.00000004 

.0070 

-.00000000 

.00000006 

.00000006 

.o")80 

-.00000000 

.00000009 

.00000009 

.0090 

.OOOO^OIC 

.00000012 

.00000002 

22.31302400 

.0100 

.00000010 

.00000017 

.00000007 

68.75311900 

.0110 

.00000020 

.00000022 

.00000002 

11.08661700 

.0120 

.00000020 

.00000029 

.00000009 

45.49851100 

.0130 

.00000030 

.00000036 

.00000007 

22.69731600 

.0140 

.000o0040 

.00000045 

.00000006 

14.43896200 

.0150 

.00000049 

.00000056 

.00000006 

12.62958800 

.0160 

.00000059 

.00000068 

.00000008 

13.82990400 

.0170 

.00000079 

.00000081 

.00000002 

2.15446370 

.0180 

.00000089 

.00000096 

.00000007 

7.92134280 

.0190 

.00000109 

.00000113 

.00000004 

3.71044550 

.0200 

.00000129 

.00000132 

.00000003 

2.23911960 

.0210 

.00000148 

.00000152 

.00000004 

2.55020530 

.0220 

.00000168 

.00000175 

.00000007 

4.08114530 

.0230 

.00000198 

.00000199 

.00000002 

.95648918 

.0240 

.00000227 

.00000226 

.00000000 

.30489807 

.0250 

.00000257 

.00000256 

.00000000 

.38535237 

.0260 

.00000286 

.00000287 

.00000001 

.49235507 

.027C 

.00000326 

.00000322 

.00000004 

1.22735300 

.0280 

.00000355 

.00000358 

.00000004 

1.07436710 

.0290 

.00000394 

.00000398 

.00000004 

1.02347600 

.^300 

.00000443 

.00000440 

.00000003 

.68238169 

.0310 

.00000482 

.00000485 

.00000003 

.67549689 

.0320 

.00000541 

.00000533 

.00000008 

1.46435570 

.0330 

.00000580 

.00000585 

.00000004 

.76658743 

.0340 

.00000639 

.00000639 

.00000000 

.03473345 

.0350 

.00000687 

.00000696 

.00000009 

1.30767190 

.03  -,Q 

.00000756 

.00000757 

.00000001 

.14157142 

.03/0 

.00000825 

.00000821 

.00000003 

.38265613 

.0380 

.00000893 

.00000889 

.00000004 

.41163880 

.0390 

.00000961 

.00000960 

.00000000 

.02896065 

.0400 

.00001029 

.00001035 

.00000007 

.64514467 

.0410 

.00001117 

.00001114 

.00000003 

.23417417 

.0420 

.00001194 

.00001197 

.00000003 

.21690074 

.0430 

.00001282 

.00001283 

.00000002 

.13668736 

.0440 

.00001369 

.00001374 

.00000005 

.34916144 

.0450 

.00001466 

.00001469 

.00000002 

.16359680 

.0460 

.00001563 

.00001568 

.00000005 

.29545440 

.0V70 

.00001670 

.00001671 

.00000000 

.03774884 

•  Oh80 

.0000^.777 

.00001778 

.00000002 

.10648549 

.0490 

.0000i893 

.00001891 

.00000002 

.12265103 
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Table  10. 

(cont. ) 

TIME 

EXACT 

APPRCX. 

ERROR 

REL.  ERROR 

.0500 

.00001999 

.00002007 

.00000008 

.42159193 

.0510 

.00002125 

.00002128 

.00000004 

.17563834 

.0520 

.00002250 

.00002254 

.00000004 

.17418449 

.0530 

.00002386 

.00002385 

.00000000 

.02488138 

.0540 

.00002521 

.00002521 

.00000000 

.00808409 

.0550 

.00002666 

.00002661 

.00000004 

.16040694 

.056C 

.00002800 

.00002807 

.00000007 

.26088540 

.0570 

.00002964 

.00002958 

.00000006 

.21516400 

.  o  ">  8  0 

.00003108 

.00003114 

.00000006 

.20077803 

.0590 

.0000T272 

.00003276 

.00000004 

.11684211 

.0600 

.00003445 

.00003442 

.00000003 

.07769276 

.0610 

.00003618 

.00003615 

.00000003 

.09519269 

.062  0 

.00003790 

.00003792 

.00000002 

.05626995 

.0630 

.00003973 

.00003976 

.00000003 

.08475895 

•  0640 

.00004164 

.00004165 

.00000001 

.02860979 

.0650 

.00004366 

.00004360 

.00000006 

.13060567 

.0660 

.00004567 

.00004561 

.00000006 

.12637703 

.0670 

.00004-767 

.00004768 

.00000000 

.01897582 

.0680 

.00004988 

.00004981 

.00000007 

.13416859 

.0690 

.00005197 

.00005200 

.00000003 

.06409186 

.0700 

.00005426 

.00005425 

.00000000 

.01633710 

.0710 

.00005655 

.00005657 

.00000002 

.02841517 

.0720 

.00005894 

.00005895 

.00000001 

.02310950 

.0730 

.00006142 

.00006140 

.00000002 

.03323103 

.0740 

.00006389 

.00006391 

.00000001 

.02234434 

.0750 

.00006646 

.00006648 

.00000003 

.03786248 

.0760 

.00006912 

.00006912 

.00000000 

.01033167 

.0770 

.00007178 

.00007184 

.00000006 

.08429641 

.0780 

.00007463 

.00007462 

.00000001 

.01830787 

.0790 

.00007748 

.00007746 

.00000001 

.01803470 

.0800 

.00008032 

.00008038 

.00000007 

.08474703 

.0810 

.00008335 

.00008337 

.00000003 

.03039534 

.0820 

.00008638 

.00008644 

.00000006 

.06947233 

.0830 

.00008961 

.00008957 

.00000004 

.04000201 

.0840 

.00009271 

.00009278 

.00000006 

.06897400 

.0850 

.00009603 

.00009606 

.00000003 

.03484085 

.  0  8  ')  C 

.00009943 

.00009941 

.00000002 

.01553851 

.08/0 

.00010282 

.00010285 

.00000002 

.02381779 

.0880 

.00010631 

.00010635 

.00000004 

.03775580 

.0890 

.00010989 

.00010994 

.00000005 

.04209515 

.0900 

.00011357 

.00011360 

.00000004 

.03329359 

.0910 

.00011723 

.00011734 

.00000012 

.09813205 

.0920 

.00012119 

.00012117 

.00000003 

.02400266 

.0930 

.00012505 

.00012507 

.00000002 

.01701779 

.0940 

.00012899 

.00012905 

.00000006 

.04744642 

.0950 

.00013312 

.00013311 

.00000000 

.00682084 

.0960 

.00013725 

.00013726 

.00000001 

.00843737 

.o->70 

.00014146 

.00014149 

.00000002 

.01670396 

.0980 

.0001V577 

.00014580 

.00000003 

.01873466 

.0990 

.0001j017 

.00015020 

.00000003 

.01714706 
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Table  11. 
Y/D  Y=1>MQDIFIED    SIMPSONS    CONVOLUTION 


TIME 

EXACT 

APPROX. 

ERROR 

REL.  ERROR 

.060 

.05992810 

.05992200 

.00000610 

.01017886 

.090 

.08^75780 

.08974500 

.00001280 

.01426060 

.120 

.11942730 

.11941100 

.00001630 

.01364847 

.150 

.14888500 

.14886500 

.00002000 

.01343319 

.180 

.17808090 

.17805700 

.00002390 

.01342087 

.210 

.20696650 

.20693800 

.00002850 

.01377034 

.240 

.23549580 

.23546300 

.00003280 

.01392806 

.270 

.26362490 

.26358500 

.00003990 

.01513514 

.300 

.29131260 

.29126800 

.00004460 

.01531001 

.330 

.31852080 

.31847000 

.00005080 

.01594872 

.360 

.34521410 

.34515700 

.00005710 

.01654046 

.390 

.37136030 

.37130200 

.00005830 

.01569904 

.420 

.39693050 

.39686900 

.00006150 

.01549390 

.450 

.42189900 

.42183400 

.00006500 

.01540653 

.480 

.44624360 

.44617500 

.00006860 

.01537277 

.510 

.46994520 

.46987500 

.00007020 

.01493791 

.540 

.49298800 

.49291700 

.00007100 

.01440197 

.570 

.51535930 

.51528800 

.00007130 

.01383501 

.600 

.53704960 

.53697800 

.00007160 

.01333210 

.630 

.55805230 

.55797700 

.00007530 

.01349336 

.660 

.57836350 

.57828600 

.00007750 

.01339988 

.690 

.59798200 

.59790200 

.00008000 

.01337833 

.720 

.61690930 

.61682900 

.00008030 

.01301650 

.750 

.63514900 

.63506900 

.00008000 

.01259547 

.780 

.65270670 

.65262800 

.00007870 

.01205748 

.810 

.66959030 

.66951100 

.00007930 

.01184306 

.840 

.68580910 

.68572800 

.00008110 

.01182545 

.870 

.70137420 

.70129400 

.00008020 

.01143470 

,900 

.71629790 

.71621700 

.00008090 

.01129418 

.930 

.73059390 

.73051200 

.00008190 

.01121006 

.960 

.74427690 

.74419500 

.00008190 

.01100397 

.990 

.75736240 

.75728000 

.00008240 

.01087986 

1.020 

.76986660 

.76978600 

.00008060 

.01046935 

1.050 

.78180640 

.78172500 

.00008140 

.01041179 

1.080 

.79319910 

.79311700 

.00008210 

.01035049 

1.110 

.80406240 

.80398200 

.00008040 

.00999922 

1.140 

.81441410 

.81433500 

.00007910 

.00971250 

1.170 

.82427220 

.82419400 

.00007820 

.00948716 

1.200 

.83365460 

.83357600 

.00007860 

.00942837 

1.230 

.84257940 

.84250300 

.00007640 

.00906739 

1.260 

.85 -.06410 

.85098700 

.00007710 

.00905925 

1.290 

.85912660 

.85905200 

.00007460 

.00868324 

1.320 

.86678400 

.86670900 

.00007500 

.00865267 

1.350 

.87405330 

.87398000 

.00007330 

.00838622 

1.380 

.88095130 

.88087600 

.00007530 

.00854758 

1.410 

.68749420 

.88741800 

.00007620 

.00858597 
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Table  11.     (cont.) 


T 

rME 

EXACT 

APPRCX. 

ERROR 

REL.  ERROR 

1 

.440 

.89369780 

.89362100 

.00007680 

.00859351 

1 

.470 

.89957750 

.89950200 

.00007550 

.00839283 

1 

.500 

.90514830 

.90507200 

.00007630 

.00842956 

1 

.530 

.91042460 

.91035000 

.00007460 

.00819398 

1 

.560 

.91542050 

.91534500 

.00007550 

.CU824758 

1 

.590 

.92014940 

.92007500 

.00007440 

.00808564 

1 

.620 

.92462430 

.92455200 

.00007230 

.00781939 

1 

,650 

.92885770 

.92878800 

.00006970 

.00750384 

1 

.680 

.93286160 

.93279100 

.00007060 

.00756811 

1 

.710 

.93664760 

.93657800 

.00006960 

.00743076 

1 

.740 

.94022670 

.94015500 

.00007170 

.00762582 

1 

.770 

.94360950 

.94353800 

.00007150 

.00757729 

1 

,800 

.94680610 

.94673500 

.00007110 

.00750946 

1 

.830 

.94982610 

.94975700 

.00006910 

.00727502 

1 

.860 

.95267890 

.95261100 

.00006790 

.00712727 

1 

.890 

.95537320 

.95530600 

.00006720 

.00703390 

1 

.920 

.95791740 

.95784900 

.00006840 

.00714049 

1 

.950 

.96031940 

.96025000 

.00006940 

.00722676 

\ 

.980 

.96258700 

.96251900 

.00006800 

.00706430 

2 

.010 

.96*72740 

.96466100 

.00006640 

.00688277 

2 

.040 

.96674730 

.96668200 

.00006530 

.00675461 

2 

.070 

.96865350 

.96859000 

.00006350 

.00655549 

2 

.100 

.97045200 

.97038900 

.00006300 

.00649182 

2 

.130 

.97214880 

.97208500 

.00006380 

.00656278 

2 

.160 

.97374940 

.97368700 

.00006240 

.00640822 

2 

.190 

.97525920 

.97519800 

.00006120 

.00627525 

2 

.220 

.97668320 

.97662000 

.00006320 

.00647088 

2 

.250 

.97802620 

.97796500 

.00006120 

.00625750 

2 

.280 

.97929260 

.97923100 

.00006160 

.00629025 

2 

.310 

.98048670 

.98042800 

.00005870 

.00598682 

2 

.340 

.98161260 

.98155400 

.00005860 

.00596977 

2 

.370 

.98267420 

.98261600 

.00005820 

.00592261 

2 

.400 

.98367490 

.98361600 

.00005890 

.00598775 

2 

.430 

.98461830 

.98455900 

.00005930 

.00602264 

2 

.460 

.98550760 

.98544800 

.000C5960 

.00604764 

2 

.490 

.98634580 

.98628900 

.00005680 

.00575863 

2 

.520 

.98713580 

.98708000 

.00005580 

.00565272 

2 

.550 

.98788040 

.98782500 

.00005540 

.00560797 

2 

.580 

.98858220 

.98852600 

.00005620 

.00568491 

2 

.610 

.98924360 

.98918600 

.00005760 

.00582263 

2 

.64  0 

.98986680 

.98980900 

.00005780 

.00583917 

2 

.670 

.99045410 

.99039500 

.00005910 

.00596696 

2 

.700 

.99100750 

.99094800 

.00005950 

.00600399 

2 

.730 

.99152900 

.99146900 

.00006000 

.00605126 

2 

.760 

.99202030 

.99196100 

.00005930 

.00597770 

2 

.790 

.99248330 

.99242500 

.00005830 

.00587415 

2 

.820 

.99291950 

.99286300 

.00005650 

.00569029 

2 

.850 

.99333040 

.99327300 

.00005740 

.00577854 

2 

.880 

.99371760 

.99366200 

.00005560 

.00559515 

2 

.910 

.99408240 

.99402700 

.00005540 

.00557298 

2 

.940 

.99442610 

.99437200 

.00005410 

.00544032 

2 

.970 

.99474980 

.99469600 

.00005380 

.00540840 
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Table  12. 
Y=1»M0DIFIED    SIMPSONS    CONVOLUTION 


time 

EXACT 

APPROX. 

ERROR 

REL.  ERROR 

.002 

.00200000 

.00190000 

.00010000 

5.00000000 

.003 

.00300000 

.00250000 

.00050000 

16.66666700 

.004 

.00400000 

.00310000 

.00090000 

22.50000000 

.005 

.00500000 

.00370000 

.00130000 

26.00000000 

.  006 

.00599990 

.00430000 

.00169990 

28.33213900 

.007 

.00699990 

.00490000 

.00209990 

29.99900000 

.008 

.00799980 

.00550000 

.00249980 

31.24828100 

.009 

.00899980 

.00610000 

.00289980 

32.22071600 

.oio 

.00999970 

.00670000 

.00329970 

32.99799000 

.on 

.C1099960 

.00730000 

.00369960 

33.63395000 

.012 

.01199950 

.00790000 

.00409950 

34.16392300 

.013 

.01299930 

.00850000 

.00449930 

34.61186400 

.014 

.01399910 

.00910000 

.00489910 

34.99582100 

.015 

.01499890 

.00970000 

.00529890 

35.32859100 

.016 

.C1599870 

.01030000 

.00569870 

35.61976900 

.017 

.01699840 

.01090000 

.00609840 

35.87631800 

.018 

.01799810 

.01150000 

.00649810 

36.10436700 

.019 

.01899770 

.01210000 

.00689770 

36.30807900 

.020 

.01999730 

.01270000 

.00729730 

36.49142600 

.021 

.02099690 

.01330000 

.00769690 

36.65731600 

.022 

.02199650 

.01390000 

.00809650 

36.80812900 

.023 

.02299590 

.01450000 

.00849590 

36.94528200 

.024 

.02399540 

.01510000 

.00889540 

37.07127200 

.025 

.02499480 

.01570000 

.00929480 

37.18693500 

.026 

.02599420 

.01630000 

.00969420 

37.29370400 

.027 

.C2699350 

.01690000 

.01009350 

37.39233500 

.028 

.02799270 

.01750000 

.01049270 

37.48370100 

.029 

.02899190 

.01810000 

.01089190 

37.56876900 

.030 

.02999100 

.01870000 

.01129100 

37.64796100 

.031 

.03099010 

.01930000 

.01169010 

37.72204700 

.032 

.03198910 

.01990000 

.01208910 

37.79131000 

.033 

.03298810 

.02050000 

.01248810 

37.85637900 

.034 

.03398690 

.02110000 

.01288690 

37.91725600 

.035 

.03496570 

.02170000 

.01328570 

37.97465800 

.036 

.03598450 

.02230000 

.01368450 

38.02887400 

.037 

.03698320 

.02290000 

.01408320 

38.07999300 

.038 

.03798170 

.02350000 

.01448170 

38.12809900 

.039 

.03898030 

.02410000 

.01488030 

36.17389800 

.040 

.03997870 

.02470000 

.01527870 

38.21710100 

.041 

.04097700 

.02530000 

.01567700 

38.25804700 

.042 

.04197530 

.02590000 

.01607530 

38.29704600 

.  n43 

.04297350 

.02650000 

.01647350 

38.33409000 

.044 

.04397170 

.02710000 

.01687170 

38.36945100 

.  145 

.04496970 

.02770000 

.01726970 

36.40296900 

.046 

.04596760 

.02830000 

.01766760 

38.43489800 

.047 

.04696540 

.02890000 

.01806540 

36. 46533800 

.048 

.04796320 

.02950000 

.01846320 

J8.49451200 

.049 

.04896080 

.03010000 

.01886080 

36.52224600 

.050 

.04995840 

.03070000 

.01925840 

3«. 54887300 

Table  12.  (cont.) 

TIME EXACT  APPRCX. ERROR REL.  ERROR 

.051  .05095580  .03130000  .01965580  38.57421500 

.052  .05195320  .03190000  .02005320  38.59858500 

.053  .05295040  .03250000  .02045040  38.62180500 

.054  .05394760  .03310000  .02084760  38.64416600 

.055  .05494460  .03370000  .02124460  38.66549200 

.056  .05594160  .03430000  .02164160  38.68605800 

.057  .05693840  .03490000  .02203840  38.70568900 

.058  .05793510  .03550000  .02243510  38.72453800 

.059  .05893160  .03610000  .02283160  ^8. 74254200 

.060  .05992810  .03670000  .02322810  38.75994700 

.061  .06092450  .03730000  .02362450  38.77668300 

.^62  .06192070  .03790000  .02402070  38.79268200 

.063  .06291680  .03850000  .02441680  38.80807700 

.^64  .06391280  .03910000  .02481280  38.82289600 

.065  .06490860  .03970000  .02520860  38.83707200 

.066  .06590440  .04030000  .02560440  38.85082000 

.067  .06690000  .04090000  .02600000  38.86397600 

.068  .06789540  .04150000  .02639540  36.87656600 

.069  .06839070  .04210000  .02679070  38. 88870300 

.070  .06938590  .04270000  .02718590  38.90040800 

.071  .07088100  .04330000  .02758100  38.91169700 

.072  .07187580  .04390000  .02797580  38.92241900 

.073  .07287060  .04450000  .02837060  38.93284800 

.074  .07386520  .04510000  .02876520  38.94283100 

.075  .07485970  .04570000  .02915970  38.95246700 

.o76  .07585400  .04630000  .02955400  38.96169000 

.077  .07684820  .04690000  .02994820  38.97059400 

.078  .07784220  .04750000  .03034220  38.97911400 

.079  .07883610  .04810000  .03073610  38.98734200 

.080  .07982980  .04870000  .03112980  38.99521200 

.081  .08082330  .04930000  .03152330  39.00273800 

.082  .08181670  .04990000  .03191670  39.01000700 

.083  .08280990  .05050000  .03230990  39.01695300 

.084  .08380300  .05110000  .03270300  39.C2366300 

.085  .08479590  .05170000  .03309590  39.03007100 

.086  .08578860  .05230000  .03348860  39.03618900 

.087  .08678120  .05290000  .03388120  39.04209700 

.088  .08777360  .05350000  .03427360  39.04773200 

.089  .08876570  .05410000  .03466570  39.05303500 

.090  .08975780  .05470000  .03505780  39.05822100 

.091  .09074960  .05530000  .03544960  39.06309200 

.092  .09174130  .05590000  .03584130  39.06779200 

.093  .09273280  .0565000U  .03623280  39.07225900 

.094  .09372420  .05710000  .03662420  39.07656700 

.095  .09471520  .05770000  .03701520  39.08052800 

.  ")96  .09570620  .05830000  .03740620  39.08440600 

.097  .09669690  .05890000  .03779690  39.068016UO 

.098  .09768750  .05950000  .03818750  39.09149100 

.099  .09867780  .06010000  .03857780  39.09471000 
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Table  13. 
C   C   SCLL   ON  OF  DY/DT  +  Y*Y  =  1 , TRAPEZO I  DAL  CONVOLUTION 


time 

EXACT 

APPROX. 

ERROR 

REL.  ERROR 

.030 

.02999100 

.02998600 

.00000500 

.01667167 

.060 

.05992810 

.05991600 

.00001210 

.02019086 

.090 

.08975780 

.08974300 

.00001480 

.01648882 

.120 

.11942730 

.11940600 

.00002130 

.01783512 

.150 

.14888500 

.14886000 

.00002500 

.01679148 

.180 

.17808090 

.17805300 

.00002790 

.01566704 

.210 

.20696650 

.20693600 

.00003050 

.01473668 

.240 

.23549580 

.23546000 

.00003580 

.01520197 

.270 

.26362490 

.26358600 

.00003890 

.01475581 

.300 

.29131260 

.29127100 

.00004160 

.01428019 

.330 

.31852080 

.31847800 

.00004280 

.01343711 

.360 

.34521410 

.34517100 

.00004310 

.01248501 

.390 

.37136030 

.37131500 

.00004530 

.01219840 

.420 

.39693050 

.39688700 

.00004350 

.01095910 

.450 

.42189900 

.42185700 

.00004200 

.00995499 

.480 

.44624360 

.44620300 

.00004060 

.00909817 

.510 

.46994520 

.46990700 

.00003820 

.00812861 

.540 

.49298800 

.49295300 

.00003500 

.00709956 

.570 

.51535930 

.51532600 

.00003330 

.00646151 

.600 

.53704960 

.53701900 

.00003060 

.00569780 

.630 

.55805230 

.55802600 

.00002630 

.00471282 

.660 

.57836350 

.57834300 

.00002050 

.00354448 

.690 

.59798200 

.59796600 

.00001600 

.00267567 

.720 

.61690930 

.61689900 

.00001030 

.00166961 

,750 

.6351490C 

.63514400 

.00000500 

.00078722 

.780 

.65270670 

.65270400 

.00000270 

.00041366 

.810 

.66959030 

.66959200 

.00000170 

.00025389 

.840 

.68580910 

.68581600 

.00000690 

.00100611 

.870 

.70137420 

.70138300 

.00000880 

.00125468 

.900 

.71629790 

.71631100 

.00001310 

.00182885 

.930 

.73059390 

.73060800 

.00001410 

.00192994 

.960 

.74427690 

.74429400 

.00001710 

.00229753 

.990 

.75736240 

.75738200 

.00001960 

.00258793 

1.020 

.76986660 

.76988900 

.00002240 

.00290960 

1.050 

.78180640 

.78183400 

.00002760 

.00353029 

1.08  0 

.79319910 

.79323100 

.00003190 

.00402169 

1.110 

.80  *06240 

.80409900 

.00003660 

.00455189 

1.140 

•  81-+41410 

.81445700 

.00004290 

.00526759 

1.170 

.82427220 

.82432000 

.00004780 

.00579906 

1.200 

.83365460 

.83370500 

.00005040 

.00604567 

1.230 

.84257940 

.84263400 

.00005460 

.00648010 

1.260 

.85106410 

.85112100 

.00005690 

.00668575 

1.290 

.85912660 

.85918600 

.00005940 

.00691400 

1.32  0 

.86678400 

.86684500 

.00006100 

.00703751 

1.350 

.87405330 

.87411800 

.00006470 

.00740229 

1.380 

.88095130 

.88101800 

.00006670 

.00757136 

1.410 

.88749420 

.88756300 

.00006880 

.00775216 

1.440 

.89369780 

.89376800 

.00007020 

.00785500 
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rable  13.  (co 

nt.) 

T  r-'r 

EXACT 

APPRCX. 

ERROR 

REL.  ERROR 

1.470 

.89957750 

.89964900 

.00007150 

.00794818 

1.500 

.90514830 

.90522300 

.00007470 

.00825279 

1.530 

.91042460 

.91050200 

.00007740 

.00850153 

1.560 

.91542050 

.91550000 

.00007950 

.00868453 

1.590 

.92014940 

.92022900 

.00007960 

.00865077 

1.620 

.92462430 

.92470700 

.00008270 

.00894417 

1.650 

.92885770 

.92894100 

.00008330 

.00896800 

1  .680 

.93286160 

.93294600 

.00008440 

.00904743 

1.710 

.93664760 

.93673300 

.00008540 

.00911762 

1.740 

.94  ")22670 

.94031100 

.00008430 

.00896592 

1.770 

.94360950 

.94369400 

.00008450 

.00895498 

1.800 

.94680610 

.94689300 

.00008690 

.00917823 

1.830 

.94952610 

.94991300 

.00008690 

.00914904 

1.860 

.95267890 

.95276500 

.00008610 

.00903767 

1.890 

.95537320 

.95546000 

.00008680 

.00908545 

1.920 

.95791740 

.95800400 

.00008660 

.00904045 

1.950 

.96031940 

.96040900 

.00008960 

.00933023 

1  .980 

.96258700 

.96267500 

.00008800 

.00914203 

2.010 

.96472740 

.96481500 

.00008760 

.00908029 

2.040 

.96674730 

.96683500 

.00008770 

.00907166 

2.070 

.96865350 

.96874300 

.00008950 

.00923963 

2.100 

.97045200 

.97054100 

.00008900 

.00917098 

2.130 

.97214880 

.97223800 

.00008920 

.00917555 

2.160 

.97374940 

.97383700 

.00008760 

.00899615 

2.190 

.97525920 

.97534600 

.00008680 

.00890020 

2.220 

.97668320 

.97677100 

.00008780 

.00898961 

2.25C 

.97802620 

.97811500 

.00008880 

.00907951 

2. 280 

.97929260 

.97938100 

.00008840 

.00902692 

2.310 

.98048670 

.98057600 

.00008930 

.00910772 

2.340 

.98161260 

.98170100 

.00008840 

.00900559 

2.370 

.98267420 

.98276100 

.00008680 

.00883304 

2.400 

.98367490 

.98376300 

.00008810 

.00895621 

2.430 

.98461830 

.98470700 

.00008870 

.00900857 

2.460 

.98550760 

.98559600 

.00008840 

.00897000 

2.490 

.98634580 

.98643300 

.00008720 

.00884071 

2.520 

.98713580 

.98722300 

.00008720 

.00883364 

2  ,550 

.98768040 

.98796800 

.00008760 

.00886747 

2.580 

.98858220 

.98867100 

.00008880 

.00898256 

2.610 

.98924360 

.98933400 

.00009040 

.00913830 

2.640 

.98986680 

.98995800 

.00009120 

.00921336 

2.670 

.99045410 

.99054500 

.00009090 

.00917761 

2.700 

.99100750 

.99109800 

.00009050 

.00913212 

2.730 

.99152900 

.99161900 

.00009000 

.00907689 

2.760 

.99202030 

.99211000 

.00008970 

.00904215 

2.790 

.99248330 

.99257400 

.00009070 

.00913869 

2.820 

.99291950 

.99301000 

.00009050 

.00911454 

2.850 

.99333040 

.99342100 

.00009060 

.00912083 

?.880 

.99371760 

.99380900 

.00009140 

.00919778 

2.910 

•  99tu8240 

.99417400 

.00009160 

.00921453 

2.940 

.99-+42610 

.99451700 

.00009090 

.00914095 

2.970 

.99474980 

.99483900 

.00008920 

.00896708 
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7  T'-'f 

EXACT 

APPROX. 

ERROR 

REL.  ERROR 

.001 

.00100000 

.00120000 

.00020000 

20.00000000 

.002 

.00200000 

.00220000 

.00020000 

10.00000000 

.003 

.00300000 

.00320000 

.00020000 

6.66666670 

.004 

.00400000 

.00420000 

.00020000 

5.00000000 

.005 

.00500000 

.00520000 

.00020000 

4.00000000 

.006 

.00599990 

.00620000 

.00020010 

3.33505560 

.007 

.0069999C 

.00720000 

.00020010 

2.85861230 

.008 

.00799980 

.00820000 

.00020020 

2.50256260 

.009 

.00899980 

.00920000 

.00020020 

2.22449390 

•  oio 

.00999970 

.01020000 

.00020030 

2.00306010 

.011 

.01099960 

.01120000 

.00020040 

1.82188440 

.012 

.01199950 

.U1220000 

.00020050 

1.67090300 

.013 

.01299930 

.01320000 

.00020070 

1.54392930 

•  014 

.01399910 

.01420000 

.00020090 

1.43509230 

.015 

.01499890 

.01520000 

.00020110 

1.34076500 

.016 

.01599870 

.01620000 

.00020130 

1.25822720 

.017 

.01699840 

.01720000 

.00020160 

1.18599400 

.018 

.01799810 

.01820000 

.00020190 

1.12178510 

.019 

.01899770 

.01920000 

.00020230 

1.06486570 

.020 

.01999730 

.02020000 

.00020270 

1.01363680 

.021 

.02099690 

.02120000 

.00020310 

.96728565 

.022 

.02199650 

.02220000 

.00020350 

.92514718 

.023 

.02299590 

.02320000 

.00020410 

.88754952 

.024 

.02399540 

.02420000 

.00020460 

.85266343 

.025 

.02499480 

.02520000 

.00020520 

.82097076 

.026 

.02599420 

.02620000 

.00020580 

.79171507 

.027 

.02699350 

.02720000 

.00020650 

.76499898 

.028 

.02799270 

.02820000 

.00020730 

.74055021 

.029 

.02899190 

.02920000 

.00020810 

.71778669 

.030 

.02999100 

.03020000 

.00020900 

.69687573 

.031 

.03099010 

.03120000 

.00020990 

.67731308 

.032 

.03198910 

.03220000 

.00021090 

.65928707 

.033 

.03298810 

.03320000 

.00021190 

.64235285 

•  034 

.03398690 

.03420000 

.00021310 

.62700629 

.035 

.03498570 

.03520000 

.00021430 

.61253598 

.o36 

.03598450 

.03620000 

.00021550 

.59886896 

.037 

.03698320 

.03720000 

.00021680 

.58621212 

.038 

.03798170 

.03820000 

.00021830 

.57475047 

.039 

.03898030 

.03920000 

.00021970 

.56361803 

.040 

.03997870 

.04020000 

.00022130 

.55354476 

.041 

.04097700 

.04120000 

.00022300 

.54420773 

.^42 

.04197530 

.04220000 

.00022470 

.53531482 

.043 

.04297350 

.04320000 

.00022650 

.52706901 

.  o44 

.04397170 

.04420000 

.00022830 

.51919757 

.n45 

.04496970 

.04520000 

.00023030 

.51212261 

.046 

.04596760 

.04620000 

.00023240 

.50557349 

.047 

.04696540 

.04720000 

.00023460 

.49951667 

.048 

.04796320 

.04820000 

.00023680 

.49371185 

.049 

.04896080 

.04920000 

.00023920 

.48855411 

•  050 

.04995840 

.05020000 

.00024160 

.48360236 
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Table  1^.  (cont.) 

TIMF EXACT  APPRCX.      ERROR  REL.  ERROR 

.051  .05095580  .05120000  .00024420  .47923887 

.052  .05195320  .05220000  .00024680  .47504292 

.053  .05295040  .u5320000  .00024960  .47138454 

.054  .05394760  .05420000  .00025240  .46786141 

.055  .05494460  .05520000  .00025540  .46483185 

.056  .05594160  .05620000  .00025840  .46191028 

.057  .05693840  .05720000  .00026160  .45944389 

.056  .05793510  .05820000  .00026490  .45723577 

.059  .05893160  .05920000  .00026840  .45544326 

.060  .05992810  .06020000  .00027190  .45371036 

.061  .06092450  .06120000  .00027550  .45219903 

.062  .06192070  .06220000  .00027930  .45106079 

.063  .06291680  .06320000  .00028320  .45011825 

.064  .06391280  .06420000  .00028720  .44936226 

.065  .06490860  .06520000  .00029140  .44893897 

.066  .06590440  .06620000  .00029560  .44852847 

.067  .06690000  .06720000  .00030000  .44843049 

.068  .06789540  .06820000  .00030460  .44863128 

.069  .06889070  .06920000  .00030930  .44897207 

.070  .06988590  .07020000  .00031410  .44944688 

.071  .07088100  .07120000  .00031900  .45005008 

.072  .07187580  .07220000  .00032420  .45105585 

.073  .07287060  .07320000  .00032940  .45203415 

.074  .07386520  .07420000  .00033480  .45325810 

.075  .07485970  .07520000  .00034030  .45458371 

.076  .07585400  .07620000  .00034600  .45613943 

.077  .07684820  .07720000  .00035180  .45778561 

.078  .07784220  .07820000  .00035780  .45964785 

.079  .07883610  .07920000  .00036390  .46159057 

.080  .07982980  .08020000  .00037020  .46373660 

.081  .08082330  .08120000  .00037670  .46607847 

.082  .08181670  .08220000  .00038330  .46848626 

.083  .08280990  .08320000  .00039010  .47107894 

.084  .08380300  .08420000  .00039700  .47373006 

.085  .08479590  .08520000  .00040410  .47655606 

.086  .0857b860  .08620000  .00041140  .47955090 

.087  .08678120  .08720000  .00041880  .48259300 

•088  .08777360  .08820000  .00042640  .48579527 

.089  .08876570  .08920000  .00043430  .48926556 

.090  .08975780  .09020000  .00044220  .49265913 

.091  .09074960  .09120000  .00045040  .49631073 

.092  .09174130  .09220000  .00045870  .49999291 

.093  .09273280  .09320000  .00046720  .50381-311 

.094  .09372420  .09420000  .00047580  .50765971 

.095  .09471520  .09520000  .00048480  .51185026 

.096  .09570620  .09620000  .00049380  .51595403 

.097  .09669690  .09720000  .00050310  .52028555 

.098  .09768750  .09820000  .00051250  .52463212 

.099  .09867780  .09920000  .00052220  .52919704 
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SUMMARY 

Investigation  of  the  graphs  of  the  first  order  equation  shows  that 
there  exists  a  sampling  interval  size  which  minimizes  the  largest  absolute 
value  of  error  for  the  fixed  computation  duration  in  Trapezoidal 
Convolution  and  the  Modified  Simpson's  Convolution.  Simpson's  1/3  rule 
shows  a  smaller  absolute  maximum  error  than  the  other  two  but  it  is 
oscillatory  and  does  not  tend  to  give  a  desired  time  interval  which 
minimizes  the  maximum  orror. 

In  determining  a  sampling  period  which  yields  a  minimum  absolute 
maximum  error  it  has  been  shown  by  Halijak  (12),  that  T  is  a  function  of 
the  solution  time,  t.  This  conclusion  was  arrived  at  using  Trapezoidal 
Convolution  and  various  differential  equation  solutions. 

At  about  T  =  .008  seconds  the  Modified  Simpson's  Convolution  starts 
and  continues  to  give  better  results  than  Trapezoidal  Convolution.  Its 
oscillation  gives  a  definite  minimum  at  T  =  .026  seconds  and  at  T  =  .052 
the  oscillation  discontinues  and  the  curve  fits  a  pattern  similar  to  that 
of  Trapezoidal  Convolution,  but  with  better  results. 

The  tendency  of  Simpson's  formula  to  give  oscillatory  results  has 
definitely  been  ameliorated  by  using  the  Modified  Simpson's  Convolution 
over  Simpson's  1/3  rule;  this  can  be  seen  from  the  graphs.  There  also 
tends  to  be  a  pattern  when  using  the  Modified  Simpson's  Convolution.  The 
error  curve  seems  to  be  fairly  symmetrical,  on  both  sides  of  the  desired  T, 
of  .026  seconds,  and  after  a  fashion,  it  tends  to  stair  step  down  to,  and 
up  from,  T  =  .026  seconds.  This  observation  could  possibly  be  useful  in 
determining  the  absolute  maximum  error  of  other  differential  equations 
when  using  the  same  solution  time. 
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ft*  first  order  equation  was  the  only  one  used  in  determining  * 
stability  analysis  of  the  absolute  Madam  error,  but  the  results  «f 
higher  order  equations  should  follow  a  stellar  pattern,    The  second 
order,  third  order  and  non-linear  differential  equations  that  were 
solved  show  only  the  dooteal  plaoe  admracy  that  oan  be  expected  for 


In  chocking  the  second  order  equation  for  accuracy,  saapllng 
of  #0005  socon03,  .001  seconds,  and  .01  seoonds  wore  U3ed 
with  ft  eolation  ttee  of  .1  eeoond  for  the  two  aaallest  sampling  intervals, 
and  1  eeoond  for  the  largest  interval.    Staaller  intervals  xjoto  not  used 
befall  of  the  larger  nuefcer  of  iterations  required  and  thus  nuah  longer 
iMIimtun   tine*    In  omparlng  the  data  between  Trapaaoidal  Convolution 
ml  the  Modified  Stepson's  Convolution  it  was  observed  that  both  nethods 
Give  vory  mam  results.    With  T  «  .01  seoonds,  the  Modified  aiapeon's 
Convolution  had  superior  accuracy  for  the  first  half  of  the  solution 
ttee  and  Trapeaoidal  Convolution  was  superior  the  last  half.    Both 
nethods  were  never  below  a  five  dooteal  plaoe  aoouraoy  and  for  the  first 
•5  seoonds  the  aoouraoy  was  consistently  six  deoteal  places,    With 
T  ~  ,001  both  nethods  showed  an  increase  in  aoouraoy  :.1th  Trapeaoidal 
Convolution  being  consistent  with  ooven  decteoi  place  aoouraoy  and  the 
IMMrm  Stamen's  Convolution  having  seven  plaoe  aoouraoy  exoept  for  the 
last  #03  seconds  which  yielded  six  place  accuracy.    With  T  *  .0005  seoonds 
the  marall  accuracy  of  both  nethods  decreased*  with  Erapesoidal  Convolution 
showing  a  superior  overall  accuracy  of  six  and  seven  dooteal  places. 


Tb»  mm  sampling  intervals,  etc.,  were  uaod  for  the  third  order 

equation  and  the  accuracy  was  approximate^  the  sane  as  that 
obtained  for  the  sooond  order  equation  esaoopt  the  ibdifiod  aiapaon's 
Convolution  lost  a  decimal  place  accuracy  for  the  first  half  of  the 
solution  tine  for  T  i  .01  seconds. 

The  sampling  intervals  used  for  the  non-linear  equation  ware  ,001 
seconds,  .01  seconds,  and  .03  seconds.  With  T  •  .001  seconds,  Trapezoidal 
Convolution  shows  superior  accuracy  by  two  deoiaal  places,  having  three 
plaoe  accuracy.  With  T  *  ,01  seconds  and  .03  seconds  the  accuracy  of 
both  methods  was  three  and  four  decimal  places,  respectively. 

Vbv  purpose  of  high  decimal  place  accuracy  and  considering  the 
•copies! ty  of  the  Modified  Simpson's  Convolution  a3  compared  to  Trapezoidal 
Convolution,  Trapezoidal  Convolution  appears  to  be  superior  for  equations 
up  to  and  including  the  third  order.  It  should  be  noted  that  round-off 
error  has  not  been  taken  into  account  in  arriving  at  the  preceding 
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APPKHUIX  B 

Solution  of  *£  +  y  ■  1 

y(0)  «  0 

bttflttM  iMMUMMftM 

0+t  *f 

s(s+1)      s      e+1  • 

A  ■     1 

B  *  -1 

y  a  s    s+i 

(77) 


Inverting  into  the  time  domain 

y(t)  »  1-e-t   (exact  solution)  (78) 

Trapesoidal  Convolution: 
fwm  equation  (77).  dividing  through  by  a 

y  +  78^  (79) 

Z  transforming  (79)  |  using  equations  (5k)   and  (55) 

Jfaltinlying  thrcu^i  by  1~a 

£(1-8)  •  §(1+3)2$  ^t^   {^}"j°t  h  1»  N  N  •••( 
Collecting  terns  in  powers  of  a 

(1  ♦  f)2£  •  (f  -  1)zZ^«  T^ 
uhioh  the  recurrence  relation  is 
O+i^a+Cf-Dja.!-^  (80) 
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Modified  Simpson's  Convolution! 

Z- transforming  (79)  using  equations  (55)  and  (57) 

^  *  1^  (12  +  ?      12^  +  (T  "  12)Z        TT^F 
Multiplying  through  by  1-z: 

AN  +  T(-|  +  ||  -  £)2gr  ♦  T(^L      g)s  =  T^ 

[xnj=<o.  I,  1,  1,  ...> 
Collecting  terms 

(1  +  ^|T)257  +  (§T  -  1)z2Sr  -  Xz2^  =  ^ 

'   f,)4.w  .M,.  ,,,... 


From  which  the  recurrence  relation  is 


**&**++•**#*&•**'*  (81) 

Next  y^  and  y2  have  to  be  determined; 
From  the  occurrence  relation 

(1  +  frfy  *£fe  •  f  (82) 

(f1-1)^  +(1  +^)y2  =  T  (83) 

Solving  (82)  and  (83)  for  y.  and  y  we  have 
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T             -ll 

T         |*£f 

V       — 

J\  - 

1  +|r     -  -2. 

3             12 

_                   _ 

I2 

1   +  t  +  ~ 
3 


ri  +% 


If  -  1 


y2  =  - 


2T 


1  +|t 


X 

12 


§T  -  1  1   +  ^T 


1   +  T  +  3T2 
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C   C   SOLUTION  OF  DY/DT  +  Y= 1 ♦ TRAPEZO I  DAL  CONVOLUTION 
DIMENSION  TIME(999)  ,Y(999) 

1  READ,T,N 

A  =  T/(l.  +  T/2.  ) 
B  =  -(T/2.-l. ) /( 1.  +  T/2.  ) 
Yd  )=0.0 
PRE=0.u 
DO  2  1=2, N 
Z=  I  -1 

Y( I  )=A+B*Y( 1-1  ) 
X=( 1.0-EXP(-T*Z)  ) 
E  =  ABS(X-Y(  I  )  ) 
IF(E-PRE)  2,2,8 
8  PRE=E 

2  CONTINUE 
PUNCH  11,T,PRE 

11  F0RMAT(F4.3»F16.8) 
GO  TO  1 
END 
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C   C   SOLUTION  OF  DY/DT  +  Y= 1 ,MODI F I  ED  SIMPSONS  CONVOLUTION 
DIMENSION  TIME(999)  ,Y(999) 

1  READ,T,N 
A=l.+(5./12. )*T 

B  =  -( (2./3. )*T-1. ) 
C=( 1./12. )*T 
Y(l  )=0.C 

Y(2)=(T+(1./2.)*T*T)/(1.+T+(1./3.)*T*T) 
PRE=0.0 
PUNCH  10 
DO  2  1=3. N 
Z  =  I-1 

Y { I )=(T+B*Y( 1-1 )+C*Y( 1-2) ) /A 
X=( 1.0-EXP(-T*Z)  ) 
E  =  ABS(X-Y( I )  ) 
IF(E-PRE)  2,2,8 
8  PRE=E 

2  CONTINUE 
PUNCH  lltT.PRE 

10  F0RMAT(//2X1HT,6X15HABS.  MAX.  ERROR,  3X//) 

11  F0RMAT(F4.3,F16.8) 
GO  TO  1 

END 
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Solution  of  SLl  +  2&L  +  y  =  1    y(o)  =  y(o)  =  0 
dtz   dt 


Laplace  transforming 


3— f(t)  =  8"f  (a) 
dtn 


-  ^   sn"P  fH 


(0) 


2-  -      -      i 

=  s  y  +  2sy  +  y  =  s 


?  -      1 

(T  +  2s  +  1  )y  =  J 


m 


7  = 


1 


1 


J> 


s( 8*428+1)       s(s+1)' 


s       s+1        (a+1)2 


For  repeated  roots 

R(S)=  m  ^$*%m  =  -2^~l&+J^  +  ...lx— 

^s>         QW   (s+s^)1,   s+gj   (s+aj)2      (S+Sj)r 


Kjr  -  R(s) 

s=-sj 

Now  when 

s  =  0,  A  =  1 
s  ■  -1 ,  C  =  -1 


K   =  ^_  4^  R(8) 
*      (r-n)!  ds(r-n) 


S=-Sj 


(2-1 ) !   s2 


-1 


s=-1 


and 

.   1 

7  = 


k      (s+1)   (s+1)2 
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Inverting  into  the  time  domain 

L  e-stfcn  m  _nj 

(s+?t)n+1 

y(t)  =  \-emt-emtt      (exact  solution)  (85) 

Trapezoidal  Convolutions 

From  equation  (84-),  dividing  through  by  s2 

y  *  §7  +  §  -  4  (86) 

sr        s 

Z-transforming  (86),  using  equations  (5*0,   (55) .  and  (56) 

^  +  iom^  +  tfz  -m  =  tz    2  lis  =  £  iflaad 

(1-z)    (1-z)2    (1-z)2  2  1-z   2   (1-zp 

Multiplying  through  by  (1-z)  and  collecting  terms 

(T+1)2y  +  {1?<4)z7&  +  (l-T)z2Zy  -  A 

n 


{\H0»  2»  1»   1»   1»  1»   1»   —  ( 
is  ^  ^ 


From  which  the  recurrence  relation 

(T+1)yn  +  (T2-2)yI>.1  +  0-T)yn-2  i  T^  (8?) 

Modified  Simpson's  Convolution » 

Z-transforming  (86)  using  equations  (55) »  (57)*  and  (58) 

Z^(T^  (Tz  +  T?  '  iP  +  TP^  -  ^-iy2z  + 

<$n  +kiy2)z2-T|(T4f-i>a3 

1-s     V12       r       12;^  +  WjT  -  T2's        2    (1.S)3 
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Multiplying  through  by  (1-z)2  and  collecting  teras,   (f^T,  f^^T) 

(1   +  Jt)£  +  (^  -  2  +  |)z^  +  (1  -  £  .  |T)z2Zjr  +  (~  +  f)z3#  =  ^ 


The  recurrence  relation  is 


(1  ♦  fa  *  (^  -  2  ♦  !>?„_,  *  0  =  #  -  f»V2  ♦  4  +  fan.)  '  Xr 


Next,  y.  and  y2  have  to  be  deterainedj 
From,  the  recurrence  relation 


(|r2  -  2)71  +  (1  +  |  -  ^)y2  =  T2 


Solving  (89)  and  (90)  for  y,  and  y0  yields 


H  i 

T2 
2 

T?,"t 
1  +  T  -  T 

1  +  *r 

I 

T2      2 
12  "  6 

?    r     I 


7 


fLn2     T 

1  ♦  2T  +  3^  -  y 


(88) 


(39) 
(90) 
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1  +|r 

T2 
2 

rrT2-  2 

I2 

il  u 

T2      T 
12      I 

frfU, 

|***|_ 

J  ^*V.fc» 


1   +  2T  +  rT2  -  ~ 
3  9 
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r       C   SOLUTION  OF  D ( 2 ) Y/DT ( 2 ) +2DY/DT+Y= 1 > TRAPEZOI DAL  CONVOLUTION 
DIMENSION  TIME(900) »Y(900) 

1  READ,T,N 
A=T+1. 
B=-(T*T-2. ) 
C=-(l.-T) 
D=T*Y/2. 

Y( 1 )=0.0 

Y(2 )=D/A 

PUNCH  10 

DO  2  1  =  3, N 

1  =  1-1 

TIME( I )=T*Z 

Y( I )=(T*T+B*Y( I-1)+C*Y( 1-2 )) /A 

X=( l.-EXP(-T*Z)-T*Z*EXP(-T*Z)  ) 

E=ABS(X-Y( I ) ) 

PCE=100.0*E/X 

2  PUNCH  lliTIMEt  I  )»X»Y(I  )»E»PCE 

10  F0RMAT(//4HTIME,6X5HEXACT,4X7HAPPR0X.,5X5HERR0R,5X10HREL.  ERROR//) 

11  FCRMAT(F5.4,F12.8»F11.8»F11.8»F13.8) 
GO  TO  1 

END 
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C   C   SOLUTION  OF  D ( 2  )  Y/DT ( 2 ) +2DY/DT+Y= 1 ,MODI F I  ED  SIMPSONS  CONVOLUTION 
DIMENSION  TIME(90C) ,Y(900) 

1  READ,T»N 

A  =  l.t-5./6.*T 

B=-(13./12.*T*T-2.+l./2.*T) 

C=-(l.-l./6.*T*T-3./2.*T) 

D=-( l./12.*T*T  +  l./6.*T  ) 

F«l.+2«*T+4./3.*T*T-l./9.*T*T*T*T 

G=T*T*(l,/2.+(2./3. )*T-(1./6.)*T*T) 

H=2.*T*T*( l.+(2./3. )*T-(l./3. )*T*T) 

Y(  1  )=0.0 

Y(2  )=G/F 

Y(3  )=H/F 

PUNCH  10 

DO  2  1=4, N 

Z=  1-1 

TIME( I )=T*Z 

Y(  I  )=(T*T+6*Y(  1-1 )+C*Y( I-2)+D*Y< 1-3)  ) /A 

X=(1.-EXP(-T*Z)-T*Z*EXP(-T*Z)  ) 

E=ABS(X-Y( I ) ) 

PCE=10u.Q*E/X 

2  PUNCH  11,TIME(  I  )  »X»Y(  I  )  .E.PCE 

10  FORMAT (//4HTIME,6X 5 H EXACT ,4X7 HAPPROX. , 5X5HERROR , 5X10HREL.  ERROR//) 

11  FCRMAT(F5.4»F12.8tF11.8.F11.8»F13.o) 
GO  TO  1 

END 
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APPENDIX  D 

Solution  of^  +  -&i+2&+y  =  1     y(0)  =  y(0)  =  y(0)  =  0 

•F       dt2  • 

Laplaco  transforming 

s3y  +  3s2y  +  3sy  +  y  ■  g  (91 ) 

1  _      _       _        i 


y  = 


s(s3+3s2+3s+1 )       s(s+1  )3 


=  A    ,       B  CD 

•       (s+1)       (s+1)2       (s+1  )3 

From  which  A  =  1 ,  B  =  -2,  C  =  -1 ,  D  «  -1 
7  "  s  "  (s+1)"  (s+1)2  "  (s+1) 3 


Inverting  into  the  time  domain 

y(t)  =  l-e^-e-H  -  ^e'H2     (exact  solution)  (92) 

Trapezoidal  Convolution: 

From  equation  (91 ) »  dividing  through  bj   s3 

7+^+"^  +  "¥  =  ^  (93) 

G     S      S"      £ 

Z-transforming  equation  (93)  using  equations  (55)  and  (56) 

4  ♦  k  i+s#  ♦  2£a_^  +  A      i  i+s^  a  _£,      2  Ua  2  _  rsO+3)2 

1-z  (1-z)2  0-z)2    2  1-z  (1-z)2    2  1-z  ^ 

Multiplying  throu^i  by  (1-z)3  and  oollecting  terms 

(1   ♦  {ftf  +  |C  -  3  -  h  +  3T2)zZy  +  (3  -  |r  +  2-  -  3T2)z2^  +  (|t-1)z33F  =  T3Xn 

{x„j=<o,  $|  jf,  1,  1,  1,  ...I 
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Frora  "which  the  recurrence  relation  is 
(1   ♦  h)yn  ♦  (f  -  3  +  |x)yn-1  *  (3-|T  «  £)^2 

+  (2T-1)yn-3  =  T3xn  W 

Modified  Simpson's  Convolution: 

Z-transformins  (93)  using  equations  (55),   (57)  t   (58)  and  (59)  yields 

(Us)"  *         "  12  12 

+  ^1-^ 

+ L«     (t2z  +  t2z2  +  a       „  ^2  +  t^3  m  £-£)7$ 

2(1-z)3  12  *  »  W 

£,.  .if     .  c    \  2 
-  TF-2Z  +  ("iFi  +  Ty2)z 

Multiplying  through  by  (1-z)3  and  collecting  terns,  the  recurrence  relation  is 

«  +  >>r, +  (^  - 1  -  3 +  |?3'yn.i +  (3  -  3T .  -^i2  +  Irt^ 

+  <fr-<+^+f)yn_3-<£  +  T+!>^  =  Xn  (95) 


7* 


fa  •&<  *#r2>,  tt3-^  ♦*£*!*&  +^2-^2+T^2^ 

T3,   T3,    ?3,   ... 

Noxt,  y.   and  y^k&ve  to  be  defceraiinad.     From  the  recurrence  relation 

(^T2  -  2T  .  3  +  — )y,   +  (1  +  —  +  —  -  *M  =  ^T3 
Solving  squatio M  (96)  arxi  (97)  fcr  y1  and  y2  yields 


(97) 


¥ 


t-3 


T3         +  T      _  T2 

1+T     +T   -TT_ 


yi  V 


^+T  +  Tfr 


1+2T 


24 


+  i   »j£ 


y2  •,= 


^T2-2T-3+  O*  1+^1     +  T£    #  ^p2 
3  I  8  * 

1+2T  £ 

4i  -2T-3+  ~T  *^3 


1+|T  +  |^T3  +  T%.|7f 


3^-^- 


1+2T 


Id  +  I       C 


^t2-2t-3+  It3   1  +  ?  +  —  -  3Ci 


* 


T5       T6 

1      T 


|   ♦  h        F  -  3T3  ff>$ 
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)F  D(3)Y/DT(3)+3D(2 )Y/DT( 2)+3DY/DT+Y= 
TIME(900) »Y(900) 


1, TRAPEZOIDAL  CGNV. 


4 
2 

10 
11 


SOLUTION  ; 

DIMENSION 

READ»T,N 

A=l.+(3./2. )*T 

6=-(T*T*T/2.-3.-3./2.*T+3.*T*T) 

C=-(3.~3./2.*T-3.*T*T+T*T*T/2. ) 

D  =  -( (3./2. )*T-1. ) 

F=T*T*T/4. 

G=(3.M.  )*T*T*T 

Y( 1  )=0.0 

Y(2)=F/A 

Y(3  )=(G+B*Y(2)  )/A 

PUNCH  10 

DO  2  1=4, N 

Z  =  I-x 

TIME( I ) =  T*Z 

Y( I )  =  (T*T*T  +  B*Y( I-1)+C*Y( I-2)+D*Y(  1-3)  ) /A 

X=1.-1.*EXP(-T*Z)-T*Z#EXP(-T*Z )-( 1./2. ) *T*T*Z*Z*EXP ( -1*1 ) 

E  =  ABS(X-Y( I )  ) 

IF(X)  3,5,3 

PCE=100.*E/X 

IF(ABS(PCE)-1000. )  4,5,5 

PUNCH  11,TIME(  I  )  ,X,Y(  I  )  ,E,PCE 

CONTINUE 

FORMAT (//4HTIME,6X5HEXACT,4X7HAPPROX.,5X5HERROR,5X10HREL.  ERROR//) 

F0RMAT(F5.4,F12.8,F11.8,F11.8.F13.b) 

GO  TO  1 

PUNCH  11,TIME(  I  ) ,X,Y(  I  )  ,E 

GO  TO  2 

END 
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C   SOLUTION  CF  D ( 3 ) Y/DT ( 3 ) +3D ( 2 ) Y/DT ( 2 ) +3DY/DT+Y= 1 ,MCD .  SIMP.  CCNV. 
DIMENSION  TIMEI900) ,Y(900) 

1  READ»T,N 
A=l.+( 5./4. )*T 

B=-( (13./4.)*T*T-( 1./2. )*T-3.+(13./24. )*T*T*T) 

C=-(3.-3.*T-(15./4.)*T*T+(3./8. )*T*T*T) 

D=-( (5./2.)*T-l.+(3./4.)*T*T+(l./8. )*T*T*T) 

F=( 1./4. )*T+( 1./4. )*T*T+< 1./24. )*T*T*T 

G=l.*9./2.*T+5./2.*T*T-3.*T*T*T*( 1 .-T/3. -l./9.*T*T-l. / 108.*T*T*T ) 

H=T*i*T*(  l./4.  +  (  l./4.)*T+(  l./16.)*T*T*T) 

C=T*T*T*(3./2.+7./2.*T-T*T-( 1 . /3. ) *T*T*T ) 

Y(  1  )=0.0 

Y(2  )=H/G 

Y(3)=0/G 

P=(T*T*T+(3./4. )*T*Y(2)*(-1.+T-(1./6.)*T*T)) 

Q=( 1./4. )*T*Y(3)*( l.-2.*T+2./3.*T*T) 

R  =  P+Q 

Y(4)=(R+B*Y(3)+C*Y( 2  )  ) /A 

Y(5)=(T*T*T+B*Y(4)+C*Y( 3)+D*Y(2 ) ) /A 

PUNCH  10 

DO  2  1=5, N 

Z  =  I-1 

TIME( I )=T*Z 

Y( I )=(T*T*T+B*Y( I-1)+C*Y( 1-2 )+D*Y( I-3)+F*Y( 1-4) ) /A 

X=l.-l.*EXP(-T*Z)-T*Z*EXP(-T*Z)-( 1./2. ) *T*T*Z*Z*EXP ( -T*Z ) 

E  =  ABS(X-Y( I )  ) 

IF(X)  3,5,3 

3  PCE=100.*E/X 
IF(ABS(PCE)-1000. )  4,5,5 

4  PUNCH  11,TIME( I ) ,X,Y( I ) ,E,PCE 

2  CONTINUE 

10  FORMAT (//4HT I ME,6X5 HEX ACT, 4X7HAPPROX.,5X5HERR0R,5X10HREL.  ERROR//) 

11  F0RMAT(F5.4,F12.8,F11.8,F11.8,F13.8) 
GO  TO  1 

5  PUNCH  11,TIME(  I  )  ,X,Y(  I  )  ,E 
GO  TO  2 

END 


APPENDIX  E 


Solution  of  ■*  +  t'   ■  1   y(0)  ■  ° 
dt 

The  known  exact  solution  is 

«t  _  a-t        o 

y(t)  =  tanh  t  =  £- C  =  1  -  *m 

ez   ♦  e-t      1  +  r 

Trapezoidal  Convolution: 

Laplace  transformlag 
1 


sy  +  L 


f 


Dividing  (99)  through  by  s 


-   1 
y  =  ~L 


/ 


1 


Z- transforming  equation  (100)  using  equations  (5*0  and  (55) 

Tz 


.   T  1+z 


fl 


(W 


Multiplying  through  by  (1-z) 


2y0-z)  +2(1+z)  z 


Tz 
1-z 


From  whioh  the  recurrence  relation  is 
T  2   T 


12   Z  2 


TX 


fVj  1°'  1*  1*  1»  *•  '•*[ 


Modified  Simpson's  Convolution: 
Z-transforming  equation  (100)  using  equations  (55)   and  (57) 

.2 


zy 


♦  £  (J-+fz.f)Zy2  +  (£-ff) 
1-z   12   3    12  *   12 


z  = 


Tz 
(1-z)' 
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(98) 
(99) 

(100) 

(101) 
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Multiplying  through  by  (1-z)  and  collecting  terms 

From  which  the  recurrence  relation  is 
»»  "  fm   +  frS  +  Kl  -  ^n-2  =  X« 


(102) 


(y}«  o.d  ♦£,  -£,,).  t.  t.  ... 


Next,  y.  and  y?  have  to  be  determined. 

These  starting  conditions  cannot  be  calculated  as  previously 
calculated,  therefore  y.  and  yo  will  be  found  by  applying  a  series 
solution  to  the  original  differential  equation,  dy/dt  +  y    ~   1,  (7). 


*  1  -  y2    7(0)  =  0 


as sums  a  solution 


y  =  c,t  +  c  tr  +  C-t^  ... 


y'   =  c1   +  2c0t  +  3c,  t^  ... 


(103) 
(105) 


Substituting  equations  (104)  and  (105)  into  equation  (103)  yields 

2 


Oj   +  2c2t  +  3o3t&  +  ...  =  1   - 


o^t  +  c  t2  +  Cot-^  ... 


a  1  -  Q*&  -  l«|^t?  -  (c22  +  l»|«|)^  (106) 
-  (26^  +  2o2o3)t5  -  (c32  ♦  2Clc5  +  20^)  t6 


Equating  like  power  of  t  in  equation  (106)  yields 

•,  =  t.  h  .  o.  c3  =  -  i,  c^  =  0.  0j  =  fa  c6  =  0.  =?  -  jjl 
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Therefore 

3   15    105 
7    m  nT  -  -Ml3-  inT)5  -  ^(nT)7  (107) 

yn  nl    3    15     105  c    7) 

Equation  (10?)  yields 

y  =  t  -  2L  +  1^  .  ^Zr7  (108) 

'      3   '5    105 

y?  =  2T  -  &T3  +  ~T5  -  Z2^T7  (1C;) 


15         10* 
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C   C   SOLUTION  OF  DY/DT  +  Y*Y  =  1 , TRAPEZO I  DAL  CONVOLUTION 
DIMENSION  TIME(900) ,Y(900) 

1  READ>T,N 
A=-l./T 
B=2.+1./(T*T) 
C=2./T 

Y(  1  )=0.0 

PUNCH  10 

DO  2  1=2. N 

Z  =  I-1 

TIME( I )=T*Z 

Y(  I  )=A  +  SQRT( B+C*Y( I -1 ) -Y ( I -1 ) *Y ( I - 1 )  ) 

X=l.-2./( l.+EXP(2.*T*Z) ) 

E=ABS(X-Y( I ) ) 

PCE=100.0*E/X 

2  PUNCH  ll.TIME(  I  )»X»Y( I  ).E»PCE 

10  FORMAT(/MHTIME»6X5HEXACT,4X7HAPPROX.,5X5HERROR,5X10HREL.    ERROR//) 

11  FCRMAT(F5.3.F12.8.F11.8.F11.8»F13.8) 
GO    TO    1 

END 


01 


C   C   SOLUTION  OF  DY/DT  +  Y*Y  =  1 ,MOD I F I  ED  SIMPSONS  CONVOLUTION 
DIMENSION  TIME1900) ♦Y(900) 

1  READ,T,N 
A=-(6./5. )*1./T 
B=(36./25.)*(l./(T*T))+12./5. 
C=l./5. 

Y(  1  )=0.0 

Y(2  )=T~T*T*T*< 1./3. ) * ( 1 .-2 . / 5 . *T*T+ ( 5 . 6  7 / 35 . ) *T*T*T*T ) 

PUNCH  10 

DO  2  1=3. N 

2*1-1 

TIME( I >  =  T*Z 

Y( I  )=A+SQRT(B+C*( 12 . * ( 1 . / T ) *Y (  I  -1  )  -8  .*Y  (  I -1 )**2  +  Y( I-2)**2  )  ) 

X  =  l.-2./( l.  +  EXP(2.*T*2)  ) 

E=A6S(X-Y( I ) ) 

PCE=100.0*E/X 

2  PUNCH  ll.TIMEl I ) »X»Y(  I  ) »E»PCE 

10  FORMAT(//4HTIME»6X5HEXACT,4X7HAPPROX.,5X5HERROR»5X10HREL.  ERROR//) 
GO  TO  1 
END 
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The  essential  purpose  of  this  report  has  been  to  develop  a  Z-transform 
notation  for  Simpson's  rule  similar  to  that  of  Trapezoidal  Convolution 
and  use  this  development  in  the  solution  of  some  differential  equations. 

Generally  speaking,  the  Z-transform  notation  must  be  dropped  when 
using  Simpson's  formula  because  of  the  finite  number  of  starting  values 
that  are  required.     This  makes  the  finite  difference  equation  ultimately 
recursive,  while  Trapezoidal  Convolution  requires  only  one  starting  value, 
making  it  instantly  recursive.     A  method  was  found  to  overcome  this  fault 
and  Z-transforms  were  effective  in  carrying  out  nuraerioal  examples  using 
the  Modified  Simpson's  Convolution. 

In  trying  to  find  an  approximation  of  Z  l/sn    in  a  similar  manner  as 
that  used  for  Trapezoidal  Convolution,  negative  results  were  arrived  at 
and  Criswell's  exact  Z  1/s11    was  used  in  the  Modified  Simpson's  Convolution 
formula.     With  this  final  touch,   the  integrator  substitution  program  was 
developed  and  described  in  Z-transform  notation. 

An  empirical  error  analysis  of  a  first  order  differential  equation 
was  made  to  determine  the  maximum  absolute  error  for  various  sampling 
intervals,  for  a  definite  solution  time.     The  error  analysis  shows  that 
there  exists  an  optimum  sampling  range  which  minimizes  the  largest  absolute 
value  of  error  using  Trapezoidal  Convolution  and  the  Modified  Simpson's 
Convolution;  that  Trapezoidal  Convolution  has  the  best  stability  with 
poorer  results;   that  the  Modified  Simpson's  Convolution  has  better 
results  and  its  stability  is  much  improved  over  Simpson's  1/3  rule. 

When  using  higher-order  Nevrton-Cotes  quadratures  such  as  the  Modified 
Simpson's  Convolution,   the  estimation  of  starting  points  plays  a  vital  role. 
The  better  the  estimate  the  higher  the  accuracy  in  numerical  analysis. 


Comparing  decimal  place  accuracy  between  the  Modified  Simpson's 
Convolution  and  Trapezoidal  Convolution  for  various  sampling  intervals 
it  was  observed  that  Trapezoidal  Convolution  yielded  overall  better 
accuracy  and  required  less  computer  time;  the  algebra  in  determining  the 
starting  conditions  for  the  Modified  Simpson's  Convolution  is  long  and 
tedious  for  the  second  and  third  order  differential  equations  solved; 
both  methods  dropped  in  accuracy  when  solving  the  non-linear  equation 
which  required  a  square  root  operation,  and  that  the  starting  conditions 
for  the  non-linear  equation  had  to  be  determined  by  using  a  Taylor  series 
or  a  differential  equation  series  instead  of  using  determinants  which 
were  used  for  the  other  differential  equations  solved. 

In  conclusion,  if  one  is  looking  for  stability  and  minimum  absolute 
maximum  error,  the  Modified  Simpson's  Convolution  may  suffice,  but  if 
one  is  looking  for  decimal  place  accuracy,  Trapezoidal  Convolution  is 
advantageous  up  to  and  including  the  third  order  differential  equation 
solved. 


